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tions, they establish both a functional law of large numbers and functional 
central limit theorem for this process. We build on that work here. For this, 
we consider a more detailed description of the system state given by K fi- 
nite, nonnegative Borel measures on the nonnegative quadrant, one for each 
job class. For each time and job class, the associated measure has a unit 
atom associated with each job of that class in the system at the coordinates 
determined by what are referred to as the residual virtual waiting time and 
residual patience time of that job. Under mild conditions, we prove a func- 
tional law of large numbers for this measure-valued state descriptor. This 
yields approximations for related processes such as the queue lengths and 
abandoning queue lengths. An interesting characteristic of these approxi- 
mations is that they depend on the deadline distributions in their entirety. 
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1. Introduction 

Here we consider a single-server queue fed by K arrival streams, each corre- 
sponding to a distinct job class. Upon arrival, each job declares its service time 
and deadline requirements. If a job doesn't enter service within the deadline 
time of its arrival, it abandons the queue before initiating service. Otherwise it 
remains in queue until it receives its full service time requirement. Nonabandon- 
ing jobs are served in a nonidling, first-in-first-out (FIFO) fashion. The queue 
is assumed to be overloaded, i.e., the offered load exceeds one. In addition, it 
is assumed that the interarrival times, service times, and deadline times are 
mutually independent sequences of independent and identically, but generally 
distributed random variables. One aim of this work is to provide approximations 
for functionals such as the queue-length process. 

Queueing systems with abandonment are observed in many applications. 
Hence, it is a natural phenomenon to study. For FIFO single-server queues with 
abandonment, the earliest analysis focused on models with exponential aban- 
donment [1], which is not an ideal modeling assumption. Thereafter, general 
abandonment distributions were considered in [3], which restricted to Marko- 
vian service and arrival processes. Shortly thereafter, stability conditions were 
determined in [2] without the Markovian restriction. In the last decade, there 
has been considerable progress with analyzing the critically loaded FIFO single- 
server queue with general abandonment [7, 14] . Here we focus on general aban- 
donment for the multiclass overloaded FIFO single-server model describe above. 

In [11], Jennings and Reed study the workload process associated with this 
queue under the assumption that the abandonment distributions are continu- 
ous. Note that the workload process is well defined since each job declares its 
deadline upon arrival. In particular, the value of the workload process is in- 
creased by a job's service time at its arrival time if and only if the value of 
the workload process immediately before the job's arrival is strictly less than 
the job's deadline time. As usual, the value of the workload decreases at rate 
one while the workload is positive. Under mild conditions, they establish both a 
functional law of large numbers and a functional central limit theorem for this 
process. 

In this paper, we build on that work. A goal is to provide a fluid approxi- 
mation for the vector-valued queue-length process. Because the abandonment 
distributions are not necessarily exponential, the queue-length vector together 
with the vector of residual interarrival times, class in service, and residual service 
time does not provide a Markovian description of the system state. Additional 
information about the residual deadline times is also needed. We track this in- 
formation using a measure- valued state descriptor. The state Z{t) of the system 
at time t consists of K finite, nonnegative Borel measures on one for each 
class, where M + denotes the nonnegative real numbers. Each measure consists 
entirely of unit atoms, one corresponding to each job of that class in the sys- 
tem. The coordinates of each atom are determine by two quantities associated 
with the job. The first coordinate is the residual virtual waiting time of that 
job, which is the amount of work in the system associated with this job and 
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the jobs that arrived ahead of it. The second coordinate is the residual virtual 
patience time. This is given by the job's deadline minus the time in system if it 
will abandon before entering service, and is given by the job's deadline plus its 
service time minus the time in system otherwise. Note that the atoms associated 
with jobs that will eventually be served are initially located in M.+ above the line 
of slope one intersecting the origin, while the atoms associated with jobs that 
won't be served are initially placed on or below this line. With this description, 
new jobs arrive and the corresponding unit atoms are added at the appropriate 
coordinates of the system state. Further, each coordinate of each atom decreases 
at rate one until the unit atom reaches one of the coordinate axes and exits the 
system. Jobs associated with atoms that hit the vertical coordinate axis exit 
due to service completion. Jobs associated with atoms that ultimately hit the 
horizontal coordinate axis exit due to abandonment (see Figure 1). 

Our choice of state descriptor is reminiscent of the one used by Gromoll, 
Robert, and Zwart [9] to analyze an overloaded processor sharing (PS) queue 
with abandonment. One distinction is that their work is for a single-class queue, 
and therefore has one coordinate rather than K . Another is that the first coor- 
dinate for the unit atom associated with a given job in the PS system is simply 
the job's residual service time. We use the waiting time to determine the first 
coordinate since it captures the order of arrivals, which is needed for FIFO. The 
evolution of the state descriptor in [9] is slightly more complicated than the one 
used here. In both cases, the patience times decrease at rate one, but in [9] the 
service times decrease at rate one over the number of jobs in the queue. Hence, 
the nice relationship with the line of slope one intersecting the origin present in 
the FIFO model is not present in the PS model. However, in both models it is 
true that hitting a coordinate axis corresponds to exiting the system with the 
vertical axis being associated with service completion and the horizontal axis 
being associated with abandonment. So this work provides an example of how 
the modeling framework developed for analyzing PS with abandonment can be 
adapted to yield an analysis in another abandonment setting. 

We develop a fluid approximation for this measure-valued state descriptor, 
which yields fluid approximations for the queue-length vector and other func- 
tionals of interest. These approximations depend on the entire abandonment dis- 
tribution of each job class. We begin by introducing an associated fluid model, 
which can be viewed as a formal law of large numbers limit of the stochastic 
system. Hence fluid model solutions are if-dimensional measure-valued func- 
tions with each coordinate taking values the space of finite, nonnegative Borel 
measures on that satisfy an appropriate fluid model equation (see (26)). We 
analyze the behavior of fluid model solutions. Through this analysis we identify 
a nonlinear mapping from the fluid workload to the fluid queue-length vector 
(see (30)). Using the nature of fluid model solutions, this mapping is refined 
to yield the approximations for the number of jobs of each class in queue that 
will and will not abandon (see (31) and (32)). In a similar spirit, we obtain 
approximations for the number of jobs of each class in system of a certain age 
or older (see (33) and (34)). In addition, we characterize the invariant states for 
this fluid model (Theorem 3.3). Next we justify interpreting fluid model solu- 
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tions and functionals derived from them as first order approximations of prelimit 
functionals in the stochastic model by proving a fluid limit theorem (Theorem 
3.2). It states that under mild assumptions the fluid scaled state descriptors for 
the stochastic system converge in distribution to measure- valued functions that 
are almost surely fluid model solutions. 

The paper is organized as follows. We conclude this section with a listing of 
our notation. The formal stochastic and fluid models are given Section 2.2. Then, 
in Section 3, we present the main results (Theorems 3.1, 3.2 and 3.3) and several 
approximations derived from our fluid model. Theorems 3.1 and 3.3 are proved 
in Section 4. In the final two sections, we provide the proof of Theorem 3.2. The 
proof of tightness is presented in Section 5, while the characterization of fluid 
limit points as fluid model solutions is presented in Section 6. For this, we prove 
a functional law of large numbers for a sequence of measure- valued processes 
that we refer to as residual deadline related processes (see Lemma 5.3). Since the 
residual deadline related processes don't include information about the service 
discipline, this result may be of more general interest for analyzing queues with 
abandonment. 

1.1. Notation 

The following notation will be used throughout the paper. Let N denote the 
set of strictly positive integers and let R denote the set of real numbers. For 
x, y € M, x V y = max(x, y), x A y = xmn(x,y), and x + — x V 0. For x € R, 
[x\ denotes the largest integer less than or equal to x and \x~\ denotes the 
smallest integer greater than or equal to x. For x G R, \\x\\ = \x\ and for iel 2 . 



The nonnegative real numbers [0, oo) will be denoted by R+. Let B\ denote 
the Borel subsets of R+ and B2 denote the Borel subsets of R+. On occasion 
we will use the notation R]j_ in place of M + . For a Borel set B £ Bi, i = 1,2, 
we denote the indicator function of the set B by 1b- For i = 1, 2, B € Bi and 
x £ Rj_, let B x denote the x-shift of the set B, which is given by 



Given x € R+ and B g B2, we adopt the following shorthand notation: 



For i = 1,2, B € Bi, and k > 0, let B K denote the K-enlargement of B, which is 
given by 



Notice that given i = 1,2 and B e B h (B K ) X C (B X ) K , but {B K ) X and {B X ) K 
are not necessarily the same set. In particular, the set resulting from shifting 
before enlarging may contain additional points. We adopt the convention that 
B x = (B X ) K , i.e., B x is the larger of the two sets. 



Nl = 



B x = {y£R\:x-y&B}. 



(1) 
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For i = 1,2, let Ci,(R+) denote the set of bounded continuous functions from 
W + to R. Given a finite, nonnegative Borel measure £ on R^, let L(£) denote 
the set of Borel measurable functions from R^ to R that are integrable with 
respect to £. For g £ L(£), we define (g, £) = J Ri gd( and adopt the shorthand 

notation ((B) = (1b, £) for B £ Bi. In addition, let \ '■ ^+ — >■ R be the identify 
map x( x ) — x fo r a H x £ R+- Given a finite, nonnegative Borel measure C on 
R_l_, if X € L(C) 5 i- e -j if (x, C) < oo, we say that £ has a finite first moment. 

For £ = 1,2, let Mj denote the set of finite, nonnegative Borel measures on 
Rj_. The zero measure in Mj is denoted by 0. For x £ R+, S x £ Mi is the 
measure that puts one unit of mass at x. Similarly, for x, y £ R+, 5( x ,y) £ M 2 is 
the measure that puts one unit of mass at (x,y). The set Mj is endowed with 
the weak topology, that is, for C", C £ M;, n £ N, we have £™ A ( as n — >■ 00 if 
and only if (g, ( n ) — > (g, () as n — > 00, for all g £ Cb(W + ). With this topology, 
Mj is a Polish space [13]. Also define 

Mf = {(Ci, . . . , Ck) : Cfc G M, for 1 < fc < A"}. 

Then Mf^ endowed with the product topology is also Polish space. Given £ G 
M^ and g £ n^L 1 L(£fc), we define the shorthand notation 

<<7,C) = «<?,Ci),.--,<<?,Ck})- 

Given £ = 1,2 and ( £ M^ , it will be handy to introduce the notation £ M, 
for the superposition measure, which is given by 

K 

C+(B) = ^2C k (B), for all Be Bi. (3) 
fe=i 

Let M^o denote the subset of Mi containing those measures that assign zero 
measure to the set {0}. Similarly, let M 2 ,o denote the subset of M 2 containing 
those measures that assign measure zero to the set 

C = R+ x {0}U{0} x R + . (4) 

For x 7 y £ R + , let 5+ £ M^o and S^ x ^ £ M 2 ,o respectively denote the measures 
that put a unit mass at the point x if x > and (x, y) if x, y > 0, and are the 
zero measure otherwise. For £ = 1,2, the set Mf is defined analogously to Mf 
except that Cfc £ M^o for 1 < fc < In addition, let Si,o denote those sets in 
B\ that do not contain zero. Similarly, let S 2j o denote those sets in £? 2 that do 
not meet the set C. 

Finally, we will use to denote convergence in distribution of random 

elements of a metric space. Following Billingsley [4] , we will use P and E respec- 
tively to denote the probability measure and expectation operator associated 
with whatever space the relevant random element is defined on. All stochas- 
tic processes used in this paper will be assumed to have paths that are right 
continuous with finite left limits (r.c.1.1.). For a Polish space S, we denote by 
D([0, 00), iS) the space of r.c.1.1. functions from [0, 00) into S, and we endow this 
space with the usual Skorohod Ji-topology (cf. [5]). There are six Polish spaces 
that will be considered in this paper: R, R + , Mi, Mf , M 2 , and Mf. 



Jennings and Puha/ Overloaded FIFO Queues with General Abandonment 



6 



2. The Stochastic and Fluid Models 
2.1. The Stochastic Model 

In this section, we define the model of the GI/GI/1 + GI queue serving K 
distinct customer classes, which will be used for the remainder of the paper. 

Initial condition and associated system dynamics The initial condition 
specifies the number Z + (0) of jobs in the queue at time zero, as well as the initial 
virtual sojourn time, initial patience time and class of each job. Assume that 
Z+ (0) is nonnegative integer- valued random variable that is finite almost surely. 
The initial virtual sojourn times, initial patience times and classes are the first, 
second, and third coordinates respectively of the first Z + (0) elements of the 
random sequence { (wj , pj , kj)}j^ C R+ x M + x {1, . . . , K}. For 1 < j < Z + (0), 
the initial job with initial virtual sojourn time Wj, initial patience time pj, and 
class kj is called job j. 

We assume that the elements of the sequence {wj}j^ arc finite, positive, and 
nondecreasing, which reflects the fact that the service discipline is first-in-first- 
out. In particular, the job with the smallest index is regarded as having arrived 
to the system before all other jobs, and therefore has the smallest virtual sojourn 
time. This is the job currently in service and Wi represents the time until its 
service is completed. The remaining jobs are waiting in the queue in order. For 
2 < j < Z + (0), Wj represents the amount of time that job j will remain in the 
system, provided that it doesn't abandon before entering service. 

We assume that the elements of the sequence {pj}jeN are finite and posi- 
tive, which reflects the fact that none of the jobs would be regarded as having 
abandoned the queue by time zero. In addition, we assume that wi < p\ , which 
reflects the fact that job 1 is presently in service and is therefore patient enough 
to stay in the queue until service completion. For 2 < j < Z + (0), we assume 
that pj < Wj if and only if Wj — Wj-\. When j is such that pj > Wj, job j 
is sufficiently patient to wait in the queue until service completion. Then, one 
regards Wj — Wj-i as the service time of job j and this service time is included in 
the virtual sojourn time of job j. When j is such that pj < Wj, job j abandons 
the queue at time pj before entering service. Thusly, a service time for job j is 
not included in the virtual sojourn time. 

Since the queue is nonidling and first- in-first-out, all Z+(0) jobs in the system 
at time zero will either abandon or be served by time W(0) = Wz + (o)- We refer 
to W(0) as the initial workload. 

A convenient way to express the initial condition is to define an initial random 
measure 2(0) e M^. For this, we will find it convenient to separate the sequence 

{(wj,pj,kj)}j=i°^ into K separate sequences {(wkj,pk,j)}f^}i \ one for each 
class. For 1 < k < K, let Zk(0) denote the number of class k initial jobs. Given 
1 < k < K, for 1 < j < Zk(0), let i(k,j) be the jth smallest index such that 
h(kj) = k and set Wkj — u>i(k,j) an d Pkj — Pi(k.j)- Then, for 1 < fc < K, let 
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.Zfc(O) G M 2 be given by 

z k (0) 

3=1 

which equals if Z k (0) = 0. Then let 

Z(0) = (Z 1 (0),...,Z K (0)). 
Our assumptions imply that Z(0) satishcs 

Furthermore, 

P ( max Z k (0)(C) = ) = 1. 

\l<k<K J 

In particular, Z(0) € M^ almost surely. 

Stochastic primitives and associated system dynamics The stochastic 
primitives consist of K exogenous arrival processes E k (-), 1 < k < K, K se- 
quences of service times {v k ,i}ieN, 1 < k < K, and if sequences of deadlines 
{rffe,i}ieN, 1 < k < K. We assume that the exogenous arrival processes, the 
sequences of service times, and the sequences of deadlines are all independent 
of one another. 

For a given 1 < k < K , the class fc arrival process E k {-) is a rate X k e (0, oo) 
renewal process. For t G [0, oo), E k (t) represents the number of class k jobs 
that arrive to the queue during the time interval (0,t]. We assume that the 
interarrival times are strictly positive and denote the sequence of interarrival 
times by {£fc,i}ieN- Class k jobs arriving after time zero are indexed by integers 
j > Z fc (0). For t e [0,oo), let 

A k (t) = Z k (Q) + E k (t). (6) 

Then class k job j arrives at time t k j = inf{t e [0,oo) : A k (t) > j}. Hence, for 
j' < j, tk,j> < t k ,j and we say that class k job j' arrives before class k job j. The 
inequality is strict for indices j > Z k (0). Moreover, for each j < Z k (0), t k ,j = 0. 

Given 1 < k < K, for each i £ N, the random variable v k ,i represents the 
service time of the (Z k (0) + i)th class k job. That is, class k job j > Z k (0) has 
service time v k j_z k io)- Assume that the random variables {v k ,i}ieN are strictly 
positive and form an independent and identically distributed sequence with finite 
positive mean \j \i k . Define the class k offered load to be p k — X k /^k- We assume 
that the queue is overloaded. In particular, we assume that p = J2k=i Pk> 

Given 1 < k < K, for each i e N, the random variable d k .i represents 
the deadline of the (Z k (0) + i)th class k job. That is, class k job j > Z k (0) 
has deadline d k j_ Zk (o)- Assume that the random variables {dfc,i}ieN are strictly 
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positive and form an independent and identically distributed sequence of random 
variables with common continuous distribution IV Assume that the mean 1/% 
is finite. Let -Ffc(-) denote the cumulative distribution function associated with 
Tfc, and denote its complement by Gk{ ) = 1 — -Ffc(')- 

As is the case for jobs in the system at time zero, jobs arriving after time zero 
are served in a first-in- first-out, nonidling fashion. A class k job j arriving to 
the system after time zero immediately enters service if the server is available. 
Otherwise, for class k job j to be served, it must wait until all other jobs currently 
in the queue exit via service completion or abandonment. If class k job j has not 
entered service before time tkj +c?fe,j-z fc (o)j class k job j abandons the queue at 
this time. Otherwise, when class k job j enters service, it is served for Vkj-z k (o) 
time units. 

It will be convenient to combine the exogenous arrival process and deadlines 
into a single measure- valued deadline process. 

Definition 2.1. For 1 < k < K, the class k deadline process is given by 



Then the deadline process is given by 

v(t) = (v 1 (t),...,v K (t)), ie[o,oo). 

Note that V k {-) G D([0, oo), Mi) for 1 < k < K and £>(•) e D([0, oo), Mf ). 

The workload process The workload process W(-) £ D([0, oo),R+) tracks 
as a function of time, the amount of time needed for the server to process all jobs 
currently in the system that will not abandon. If no additional jobs were to arrive 
after time t, the system would be empty W(t) time units in the future. This 
quantity also records the amount of time that a newly arriving job would have 
to wait before being served, if the job were sufficiently patient. In particular, as 
in [11], W(-) almost surely satisfies, for all t € [0, oo), 



Here -B(-) denotes the busy time process. Since the server works at rate one 
during any busy period, the amount of work served by time t is equal to B(t). 
Occasionally, it will be convenient to refer to the idle time process, which is given 
by I(t) = t — B(t) for t G [0, oo). For a given time t, the first and second terms in 
(7) add up all of the work that enters the system by time t and doesn't abandon 
before entering service. Note that because of the indicator in the integrand, a 
particular job's service time is added to the workload if and only if that job will 



E k (t) 



Vk{t) = §d ^ te [°'°°)- 



i=l 
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not abandon before it enters service. Fluid and diffusion limit results for this 
process were proved in [11], where it was referred to as the virtual waiting time 
process. We will leverage that fluid limit result to carry out the analysis here. 

Evolution of the virtual sojourn and patience times For 1 < k < K, 

let 

u>k,j, 1 < j < Zfc(O), 
W(t kj ), j>Z k (0), 

Pkj, l<J<Z k (0), 
d k ,j-z k (o) + v kj-z k (Q)^{d k . J _ Zk(0) >w(t kij -)}, 3 > Z k {0). 

Note that for each 1 < k < K and j > Z k (0), 

W(t kJ ) = W{t k ,j-) + «fc J j-Z»(0)l{d»j_* fc( o)>W r (tM-)}- 

Hence, if class k job j is served, w k j includes the job's service time in addition 
to the time spent waiting for service to begin. So, if class k job j is served, it 
stays in the system w k j time units. Therefore, w k j is referred to as the virtual 
sojourn time of class k job j. Further, p k j represents the initial patience of class 
k job j. If a class k job j arriving after time zero will enter service before time 
t k ,j + d k j_ Zki0) , the job's patience time is taken to be d k j^ Zk (o) + v k,j-z k (o) to 
account for the fact that it will not abandon d k j_z k (o) time units after arrival. 
Instead, it will stay until time t k j + w k j when it receives its full service time 
requirement. Otherwise, if class k job j won't enter service before the deadline 
expires, the job will abandon at time t k j + d k j_z k (o), an d p k ,j — d k j-z k (o)- 
Then for each 1 < k < K and j £ N, the sojourn time s k j of class k job j is 
given by 

s k ,j = mm(w k ,j,p k j). 

This quantity indicates precisely how long class k job j will reside in the system. 
For all t e [0, oo), 1 < k < K, and 1 < j < A k (t), define 

w k,j(t) = Ofcj - (t - t kJ )) + , (9) 
Pkj(t) = (Pk,J - (t ~ t kj )) + . (10) 

For l<k<K,jeN and t > t k j, w k j(t) and p k j(t) respectively represent 
the residual virtual sojourn and residual patience time of class k job j at time t. 
Then the residual sojourn time s k j(t) for class k job j at time t >t k j is given 
by 

s k ,j(t) = mm(w kt j(t),p ktj (t)). 

Measure- valued state descriptor For 1 < k < K, define the class k state 
descriptor by 

A k (t) 

= zZ S t k AtM))> ie [°'°°)- (") 



w kj = 

Pk, 3 = 
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The state descriptor is defined as 

Z(t) = (Z 1 (t),...,Z K (t)), ie[0,oc). (12) 

For each l<k<K, Z k {-) € D([0, oo), M 2 ), and Z(-) € D([0, oo), Mf ). 

Figure 1 depicts one component of a hypothetical system state at a fixed 
time. The points in the figure correspond to unit atoms of the measure. All 
points move to the left and down at rate one. The dotted diagonal line p = w 
separates the points in the figure into two groups. The jobs associated with 
points on or below the line will eventually abandon; the points above the line 
represent jobs that will be served. Once a point reaches one of the coordinate 
axes it immediately leaves the system and so is no longer included in the system 
state. Among all of the components of the system state, there is a unique point 
that has the smallest virtual sojourn time coordinate and lies above the diagonal. 
This point corresponds to the job currently in service. 

Notice that in Figure 1 there are three sets of points that are aligned verti- 
cally. In each set, at most one of these points is above the dotted line and the 
corresponding job will be served. For each virtual sojourn time assumed by some 
job that is in the system at this fixed time, there is exactly one of these jobs 
for which the location of the corresponding point in the component of the state 
descriptor associated with that job's class lies above the diagonal. This job will 
be served and it arrived before any other job with the same virtual sojourn time. 
These later arriving jobs aren't patient enough to remain in queue until service 
begins. Instead each will abandon. Therefore the corresponding points in the 
components of the state descriptor associated with those jobs' classes lie on or 
below the diagonal. Although their sequence of arrivals relative to one another 
is not captured by the state descriptor, the relative patience values reveal the 
order in which they will abandon. 

The state descriptor satisfies the following system of dynamic equations. For 
each 1 < k < K and for all B e B2.0 and t £ [0, 00), 

A k (t) 

Z k (t)(B)= J2 ^("M-ft^ ( 13 ) 
i=i 

To see this, simply note that for each B £ #2,0, x € R+, and w,p > 0, 

((10 - x) + , (p - x) + ) e B ^ (w-x,p-x)eB (w,p) e B x . 

The dynamic equation (13) is equivalent to 

A k (t) 

Z k (t)(B) = Z k (0)(B t ) + ]T lB t . tk Jw kij ,p k>j ). (14) 

j = l + Z fe (0) 

Given B £ 82,0 and x, y E M+, notice that 

(B x )y — B x -\.y. 
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Abandonment ^^^^^^^^^^ 

Virtual Sojourn Time 

Fig 1 . One coordinate of a hypothetical system state at a fixed time with depicted transitions. 



This together with (13) implies that, for each 1 < k < K and for all B 6 62,0 
and h, t € [0, 00), 

A h (t+h) 

Z k (t + h)(B) = Z k (t)(B h ) + ^B t+h ^ kj (w kj ,p kj ). (15) 

3 =A k (t) + l 

2.2. The Fluid Model 

In this section, we define the fluid model associated with this GI/GI/1+GI queue 
with K distinct job classes. The primitive data for this fluid model consists of 
the vector A of arrival rates, the vector \i of service rates, and the vector T 
of deadline distributions. The triple (A,^,T) is referred to as supercritical data 
since p > 1. We begin by summarizing the results in [11] that suggest a workload 
fluid model. Then, we develop a full measure-valued fluid model. For this, we 
first need to define the scaling regimes that yield the desired limiting dynamics. 

The Sequence of Time Accelerated Systems We consider a sequence of 
systems indexed by n S N in which the arrival rates and mean service times in 
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the nth system are sped up by a factor of n. We use a superscript n to denote 
all processes and parameters associated with the nth system. 

Specifically, the interarrival times j}ieN for class k in the nth system are 
given by ^ = £,k,i/ n f° r i € N. Then denotes the class fc exogenous 

arrival processes in the nth system. Hence, for 1 < k < K, E£(t) = Ekint) for 
t G [0,oo) and A£ = n\}~- Similarly, the class k service times {vki\ien in the 
nth system are given by = i>fe,j/n for i G N. Then /u£ = n/Zfc and p^ = pf.. 
Hence, p™ = p. 

The deadline sequence is unsealed. In particular, for each 1 < k < K and 
n, i G N, djjj = cJ/c^. Hence, we omit the superscript when referring to the class 
k deadlines, distribution T^, mean I/7/5, and cumulative distribution function 
Fk or its complement Then the class fc deadline process for the nth system 
is given by 

*??(*)= £ ^[0,oo), 

andD»(.) = (Df(.),..,,D» K (.)). 

For each n G N, there is an initial condition Z n (0) and fc")} Jg N 
satisfying the conditions indicated above. The initial conditions may vary with 
n. Then for n G N, 1 < k < K, and t € [0,oo), = Z£(0) + If 

we suppose for simplicity that Z + (0) = in the unsealed system, then the job 
arrival times {i£ jjjeN for class fc jobs in the nth system would be given by 

,„ = /o, l<i<^(0), 
W |^EW j j>^(0). 

The workload busy time B n (-), and idle time /"(•) processes for the 

nth system satisfy equations analogous to (7) and (8). The residual virtual 
sojourn w^j(-) and residual patience Pkj(') times, 1 < k < K and j £ N, for 
the nth system are defined as in (9) and (10). The class k state descriptor Z%{-) 
and the state descriptor Z n {-) for the nth system are defined as in (11) and 
(12). Analogs of (13) and (15) hold for the nth system as well. 

The Workload Fluid Model For the time accelerated scaling regime, the 
authors of [11] identify what one might refer to as a workload fluid model. 
Namely, define a workload fluid model solution to be a function w : [0, 00) — > K + 
satisfying 

w{t) =«j(0)+ Vft / G k (w(s))ds - 1, te[0,oo). (16) 
fc=l J ° 

Equation (16) can be interpreted as a fluid analog of (7) and (8). Notice that 
any solution to (16) is necessarily Lipschitz continuous with Lipschitz constant 
p—1, and therefore is almost everywhere differentiable. In [11], it is asserted that 
for each wq G a unique workload fluid model solution w(-) with w(0) = Wq 
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exists. They further show that workload fluid model solutions satisfy a nice 
monotonicity property [11, Theorem 1]. To describe this, let 



wi = sup < u € 



w u = sup < u £ 



5> fc G fc (u)<l] 
fe=i J 

Y,PkGk{u) < ll 
fe=i J 



Note that by continuity of Gk(-) for each fc and the fact that p > 1, < Wi < 
w u < oo. They show that for a workload fluid model solution w(-) such that 
w (0) [w;,ui u ], w(-) is strictly monotone and 

lim4)=( W1 ' iMOXu*, (17) 
if w(0) > w„. 

Otherwise, if iu(0) € [w;, iu u ], then w{t) = w(0) for all t £ [0, oo). 

Workload fluid model solutions also satisfy a useful relative ordering property. 
In particular, given two workload fluid model solutions u>i(-) and W2{-) such that 
wi(0) < w 2 (0) it follows that for all t G [0,oo), 

wi(t) < w 2 {t). (18) 

To see this note that if u>i(0) < w u and wi < W2(0) then (18) follows immediately 
from the monotoncity property of workload fluid model solutions. Otherwise, 
W2(0) < wi or wi(0) > w u . If wi(0) > w u , then by continuity and (17) there 
exists t £ [0,oo) such that 102 (*) = wi(0). For s £ [0, 00), set w(s) — w%{t + s). 
Then, it is straightforward to verify that w(-) is a workload fluid model solution 
with w(Q) = wi(t). Hence, by uniqueness, w(-) = Wi(-+t). Then for all solution 
with w(Q) = w 2 (t). Hence, by uniqueness, w(-) = w 2 {- +t) = Wx(-). Then for all 
s £ [0,oo), wi(s) = w 2 (s + t) < w 2 (s). An analogous argument demonstrates 
(18) when ^(O) < wt. 

In [11], the authors justify interpreting the fluid model as a first order approx- 
imation of the stochastic system by showing that if W n (0) — > wo almost surely 
as n —> 00, where wo is a finite nonnegative constant, then (W n , I n ) {w(-), 0), 
as n — > 00, where w(-) is the unique workload fluid model solution such that 
w(0) — wo [11, Theorem 1]. Here, we wish to cite a slightly modified version of 
[11, Theorem 1], which we state next. Before this, we make note of a represen- 
tation for the workload process developed in [11] that will be of use here. By 
[11, (19) and arguments presented in the proof of Theorem 1], for all n £ N and 
< s < t < 00, 

W n (t) = W n (s)+I n (t)-r(s)~(t-s) (19) 

+X n (t)-X n (s) + Y / Pk G k (W n (u))du, 
k=i Js 
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where {X n (-)} n ^ C D([0, oo),R) is such that as n — > oo, 

X n {-) ^ 0. (20) 

Theorem 2.2. If W n (0) =^ Wo as n — >• oo, tten bji-) oo, 

(W™,/™) (W*(-),0), (21) 

where W*(0) is egual m distribution to Wo- Furthermore, W*{-) is almost surely 
a workload fluid model solution. 

Summary of proof. Let I : [0, oo) — > [0, oo) be given by i(i) = t for all t G [0, oo). 
Since the limit in (20) is deterministic, it follows that as n — > oo, 

W n (0)-i(-)+X n {-)^ W -i(-). 

Using tools developed in [14], the authors further show that for all n G N, 

(rO, /"(■)) = (Vd ^Ag) (V^"(0) - + X"(-)) , (22) 

where (ipg^g) '■ D([0,oo),R) — > D([0, oo), R?_) is a Lipschitz continuous func- 
tion in the topology of uniform convergence on compact sets. Then, by the 
continuous mapping theorem, it follows that as n — > oo, 

(W n (-),I n (-))^(<p G ,i>G) (Wo -t(-)). (23) 

The authors of [11] go on to show that for wo G M+ ((Pg,iPg) {wo — <<(•)) = 
(«;(•), 0), where w(-) is the unique workload fluid model solution such that 
w(0)=w . □ 

The Measure- Valued Fluid Model To develop a measure-valued fluid 
model, we add spatial scaling to the sequence of time accelerated systems. In 
particular in order to obtain a fluid scaled system, we divide space in the nth 
time accelerated system by a factor of n. Then, for n G N, 1 < k < K, and 
t G [0, oo), 

m) = ^, m) = ? m , ^d zm = ^- 

n n n 

The nth fluid scaled state descriptor in the sequence is given by 

We are interested the the limiting behavior as n approaches infinity. 

One can speculate as to how the dynamics of the system will behave in the 
limit. Class k fluid should arrive to the system at rate Xk and be distributed over 
the positive quadrant as it enters. Class k fluid arriving at time t should have 
virtual sojourn time w(t), where w(-) is an appropriate fluid workload solution, 
and patience time distributed according to Tk ■ Then the residual virtual sojourn 
and residual patience times should each decrease at rate one until at least one 
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is zero, at which time the fluid would exit from the system. In particular, at 
time t + h, fluid that arrived to (w,p) at time t should be at (w — h,p — h) 
if (w — h) A (p — h) > and should have exited the system otherwise. Such 
departures would be associated with abandonment if p — h < A (w — h) and 
service completion if w — h < and w — h < p — h. We wish to capture these 
dynamics by defining an appropriate fluid model. 

The initial measure will need to satisfy various natural conditions. To define 
these, first recall the definition of C given in (4). Then, given d G Mf , recall 
that i?+ G M2 denotes the superposition measure (see (3)). Let 

wd = sup{x G M + : 00) x K + ) > 0}. 

Define I C Mf to be the collection of -d G such that 

(LI) d+{C x ) = for all x G R%, 

(1.2) Wtf < 00, 

(1.3) max Ghiws — e) > for all e > 0. 

l<k<K 

We refer to the collection C = {C x : x e ]R^_} as the corner sets. Then 
condition (1.1) is that each coordinate of the initial measure doesn't charge 
corner sets, which is equivalent to the condition that each coordinate of the 
initial measure doesn't charge vertical or horizontal lines. Given $ g that 
satisfies (1.1) and e > there exists k > such that 

max ^ sup # fc (C£) < e. (24) 

l<k<K xeR2+ 

To see this, given x G let Pi(x) = Xi for i = 1, 2. Then, given 1/ G M2 and 
i = 1,2, let 7Ti (i/) G Mi be the projection measure such that for all / G C^M^) 

(fMv)) = (f°Pi,v)- (25) 

Then, if $ G satisfies (LI), it follows that for i = 1,2, ^(i?.).) doesn't charge 
points. Hence, for every e > there exists k > such that 

max sup 7ri(tf + )([(y - k) + , y + «]) < -, 
(cf. [8, Lemma A.l]). Since for each 1 < k < K and x G 

MC*) < < 7r 1 ( 1 ?+)([(a ;i - kJ+.u + «]) + 7T 2 ( 1 ? + )([(X 2 - K )+,X 2 + k]), 

(24) follows. 

Condition (1.2) dictates that the fluid analog of the initial workload is finite 
and condition (1.3) requires that the initial fluid workload does not exceed the 
maximal deadline. In particular, let 

dmax = max supjx G M+ : Gk(x) > 0}. 

l<k<K 
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Then, by (1.2), w$ < d max if <i max = oo. Further, by (1.3), w» < d max if d max < 
oo. This prevents there from being an initial period of time during which all 
arriving fluid eventually abandons, i.e., an initial period of time during which 
w(-) decreases at rate one. The reader will see that it is needed mathematically 
when d max < oo to prevent fluid from building up on a moving vertical line 
during (0,w$ - d max ]. 

A fluid model solution for the initial measure G I is a function £ : [0, oo) — > 
Mf such that £(0) = and for each 1 < k < K, B G B 2 , and t G [0,oo), Cfc 



where w(-) denotes the unique workload fluid model solution with w(0) = w$. 
Notice that (26) is a fluid analog of (13). 

3. Main Results and Approximation Formulas 

Here we state the two main theorems proved in this paper (Theorems 3.1 and 
3.2), which together validate approximating various performance processes via 
fluid model solutions. Then, we derive some specific approximation formulas. 
Lastly, we identify the set of invariant states for the fluid model. 

Theorem 3.1. Let i9 G I. Then there exists a unique fluid model solution £(•) 
for the data (A, fJ.,T) such that £(0) = i9. In addition, 

1. w^( t ) = w(t) for all t G [0,oo), where w(-) is the unique workload fluid 
model solution such that w(0) = w$; 

2. C+(t)(C x ) = for all xeR 2 + and t G [0,oo); 

3. C(') is continuous. 

The proof of this theorem is given in Section 4. Property 1 is that the right 
edge of the support of the fluid model solution superposition measure is equal 
to the workload fluid model solution for all time, which is true by definition 
at time zero. It implies that (1.2) holds for all time. Then by the monotoncity 
properties of workload fluid model solutions, (1.3) holds for all time as well. The 
proof of Property 1 is fairly straight forward, as the reader will see in Section 4. 
Property 2 is that the fluid model solution superposition measure doesn't charge 
corner sets for all time, i.e., that (1.1) holds for all time, and its proof is more 
involved. For this, we first prove Lemma 4.1, from which Property 2 follows by 
letting e decrease to zero. Lemma 4.1 is then used together with some additional 
arguments to verify Property 3, continuity. 

The next result justifies regarding fluid model solutions as first order ap- 
proximations for the measure-valued state descriptor of the original stochastic 
system. 

Theorem 3.2. Suppose that Zq — (Zq 1 , . . . , Zq k ) is a random measure in 
with superposition measure Zq + such that 



satisfies 




(26) 



(A.l) P(Z*GI) = 1, 
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(A.2) E[<pi+p 2 ,Z * + >] <oc, 
(A. 3) E [Z* + (R2J] < oo. 

Let Wq — wz* ■ Also suppose that, as n — > oo, 

(Z"(0), (pi,Z n (0)) , (p 2 ,Z n (0)) , W n (0)) (27) 
^(Z*,( Pl ,Z*),(p 2 ,Z*),W^. 

Then, as n oo, 

Z n (-)^Z*(-), 

where Z*(0) is equal in distribution to Zq. Furthermore, Z*(-) is almost surely 
a fluid model solution. 

This result is proved via verifying tightness and then proving that any limit 
point is almost surely a fluid model solution, which are done in Sections 5 and 
6 respectively. 

Next, we use the fluid model to derive various approximation formulas. These 
include approximation formulas that demonstrate the nonlinear nature of state 
space collapse that takes place for this model. That is to say, one can recover 
various K dimensional quantities such as the fluid approximation for the queue- 
length vector from the fluid approximation for the workload process. However, 
this mapping depends on the entirety of the deadline distributions and is there- 
fore nonlinear in its dependence on the workload process. See (30), (31), and 
(32) below. In addition, we are able to approximate various age related processes 
such as the number of jobs in the system that are within a specific age range. 

To describe these, we will need to introduce the following functional. Given 
i? € I, let ((■) denote the unique fluid model solution such that C(0) — Let 
w(-) denote the unique workload fluid model solution such that w(0) = w$. 
Fluid that arrives at time s > has fluid workload w(s). Some of this fluid 
remains in the system at time t > s only if 

w(s) - (t - s) > 0. 

Recall that w(-) is continuous. Furthermore, it is strictly decreasing if w(0) > w u 
and bounded above by w u otherwise. Additionally, since p > 1, w u < d max . 
Finally, by (1.3), w(0) < d max . Hence, w(s) < d max for all s > 0. Therefore, 
w'(s) — J2k=i PkGk(w(s)) — 1 > — 1 for all s > 0. Then, as a function of 
s £ [0, oo), w(s) + s is continuous and strictly increasing. Furthermore, for 
s £ [0,t], the values range from w(0) to w(t) + t > t at time t. Hence, inf{0 < 
s < t : w(s) + s > t] is well defined, finite, and can be interpreted as the time 
at which fluid departing at time t via service completion arrived to the system. 
For t € [0,00), let 

r(t) = inf{0 < s <t: w(s) + s > t}. (28) 

Then r(t) = for all t E [0,w(0)], so that w(r(t)) + r(t) > t for all t G [0,to(0)). 
Further, for t > w(0), 

ui(r(t)) + r(t) = t. (29) 
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For t > w(0), all fluid in the system at time t arrived during the time interval 
(r(t),t\. 

Nonlinear State Space Collapse Given let ((■) denote the unique 

fluid model solution with ((0) = $ and u>(-) denote the unique workload fluid 
model solution with w(0) — w$. 

Queue-Length Vector Fluid Approximation. For 1 < k < K and t e [0,oo), 

z frt-rmfRM-l a(0)((R + )i) + Ak/ « Gfc(fl)dfl ' <<w(0) ' 

[Xk J Gk(a)da, t>w(0). 

We have chosen a as the variable of integration to suggest the interpretation age, 
or time in system. This provides a simple interpretation of this formula. At time 
t < w(0) none of the fluid that entered the system after time zero has been fully 
processed. So the integral term only needs to address departures resulting from 
expiring deadlines. Then, fluid arriving in (0, t] remains in the system at time t 
if and only if the initial deadline exceeds the current age. Hence the simple form 
of the integral. At time t > w(0) all fluid initially in the system has departed, 
and u>(r(i)) can be interpreted as the age of the fluid departing the system via 
service completion at time t. So, for t > r(i), u>(-r(t)) can be thought of as the 
age of the oldest fluid in the system at time t. 

Nonabandoning Jobs Queue-Length Vector Fluid Approximation. Let U — 
{(w,p) £ : w < p}. Then the mass present in U at time t is associated with 
fluid that doesn't abandon. For 1 < k < K and t e [0, oo), 

n k {t) = Ck{t){U) = i t „, (31 
[A fc J r(t) G k (w(v))dv, t > w(0). 

Here the variable of integration should be regarded as time, and the formula has 
the following interpretation. Fluid that has not departed the system by time t 
will not abandon prior to service completion if the initial deadline exceeds the 
fluid workload at the time of arrival. 

For t > w(0) 7 we can rewrite this result in an alternative form that high- 
lights the independence of the service time distributions to the deadline and 
interarrival time distributions. For 1 < k < K and t > u>(0), let 

w k (t) = Pk G k (w(v))dv. 

Jr(t) 

Then Wk(t) denotes the amount of fluid workload in the system at time t due 
to class k fluid. For each 1 < k < K, we have for all t > w(0), 

n k (t) = PkW k (t). 



In this form, the formula is reminiscent of linear state space collapse, but the 
dependence of Wk{-) on w(-) is nonlinear. 
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Abandoning Jobs Queue-Length Vector Fluid Approximation. Let L = {(w,p) £ 
Ki : p < w}. Then the mass present in L at time t is associated with fluid that 
abandons. For 1 < k < K and t € [0, oo), let a k {t) — C,k{t){L). Then, for 
1 < k < K and t £ [0,oo), 

M = lCk(0)(Lt) + X k S*(G k (t-v)-G k (w(v)))dv, t<w(0), 
Qk[) \X k f* (t) (G k (t-v)-G k (w(v)))dv, t>w(0). 

Note that it is easy to verify that Zk(t) = a k (t) + nk(t) for 1 < k < K and 
t £ [0,oo). 

Age Related Fluid Approximations Given d £ I, let CO) denote the 
unique fluid model solution with £(0) ~ $ and let w(-) denote the unique work- 
load fluid model solution with w(0) = w$. Fluid in the system at time t > 
that arrived at time s £ (0, t] is t — s units old at time t. Hence, for fluid in the 
system at time t > that arrives in (0,t\ to be of age at least u, it must have 
arrived by time t — u. Then s £ (0, t — it]. This fluid has residual offered waiting 
time (w(s) — (t — s)) + . Since it is in the system at time t, w(s) — (t — s) > 0, 
i.e., w(s) + s > t so that ,s £ (r(t),t — u]. For t > w(0), this is an empty 
time interval if w(r(i)) < u. Otherwise, for t > w(0) and s £ (r(t),t — u], 
w(s) -(t- s) £ (0, w(t -u)- u}. For t > w(0) and < u < t, let 

H(t, u) = [0, (w(t - it) - u) + ] x R+. 

This is the line {0} x R + if u > w(T(t)) and is a vertical stripe otherwise. Then, 
for 1 < k < K , t > w(0), and < u < t, we can interpret 

z k (t,u) = ( k (t){H{t,u)), 

to be the amount of class k fluid in the system at time t of age at least u. 

If t < w(0), some initial fluid may remain in the system at time t. We regard 
that fluid as being t units old. At time t, the residual offered waiting time of 
such fluid lies in [0,u>(0) — t\. If we consider a time < u < t < w(0) and ask 
for the total amount of fluid of age at least u, this would also include fluid that 
arrived after time zero and by time t — u. By (16), w(t — u) — u > w(0) — t. 
Hence the definition of H(t,u) and interpretation of z k (t,u) naturally extend 
to all < u < t < w(0). Then, for 1 < k < K we have the following: for 
< u < t < w(0), 

z k (t, u) = ( k (0)(H(t, u) t ) + X k I G k (v)dv, (33) 

J U 

and for t > w(0), 

(x kI :^G k{v) dv, 0<u< W (r m 

[0, w(T(t)) <U<t. V ' 

One can obtain similar approximations for the abandoning and nonabandoning 
queue length of a certain age or older by computing the measure of an appro- 
priately chosen set. 
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Invariant States An invariant state is a measure 9 € I such that the unique 
fluid model solution £(•) with initial measure 9 satisfies = 9 for all t E [0, oo). 
Here we identify the collection of invariant states. 

For this, we begin by recalling some measure theoretic background. Details 
can be found in [6, Chapter 1]. Let 

TZ = { [a, b) x [c, d) C R+ : a, c G R+, a < b < oo and c < d < oo} . 

Then TZ is an elementary family (a collection of sets that contains the emptyset, 
is closed under pairwise intersection, and such that complements of members 
of the collection can be written as finite unions of members of the collection) . 
Further note that cr(TZ) = Bi. Let 

TV = {JiLxRi ■ m G N and R t e K} . 

Then TV is an algebra (a collection of sets that contains the emptyset and is 
closed under pairwise union and relative complementation). In fact, TV is the 
algebra generated by 7Z. It is easy to see that any finite pre-measure (additive 
R+ valued function that assigns value zero to the emptyset) on 7Z extends to 
a finite pre-measure on TV . Then, by the Caratheodory extension theorem [6, 
Theorem 1.14], any finite pre-measure defined on TV uniquely extends to a finite 
Borel measure on R^_. Thus, in order to uniquely specify a finite Borel measure 
on R^_, it suffices to specify a finite pre-measure on TZ. 

Given Wi < w < w u , let 9 W G I be the unique measure that for 1 < k < K 
satisfies 

6%([w,ao) x R + ) = 0, 
and for < a < b < w and 0<c<d<ooe R + , 

/•w—a 

6%([a,b)x[c,d))=\ k / r k ([c + u,d + u))du. 

J w — b 

It is easy to verify that these relationships determine a finite pre-measure on 
TZ. Indeed, since [w, oo) x R + has measure zero, it suffices to specify each 9% on 
sets in TZ that don't meet [w,oo) x R + . Define 

J = {9 W : wi < w < w u }. 

Theorem 3.3. The set of invariant states is given by 3. 

The proof of Theorem 3.3 is given in Section 4. Note that for wi < w < w u and 
1 < k < K , the fluid approximations for queue length z™ and nonabandoning 
queue length take the form 

[>W 

z% = 9%(m 2 + ) = \ k J G k (u)du, 
nt ee 9l»{U) = \ k wG k (w). 



Jennings and Puha/ Overloaded FIFO Queues with General Abandonment 



21 



4. Properties of Fluid Model Solutions 

In this section, we prove Theorems 3.1 and 3.3. The proof of existence and 
uniqueness for Theorem 3.1 is relatively straight forward. Indeed, since the right 
hand side of (26) only depends on (A, fi, T) and (26) can be regarded as a 
definition, provided that the integral term is a well defined function taking 
values in M^. In this regard, note that £ M2 and S^r s -\ xTi e M2 for all 
s £ [0,oo) and 1 < k < K, Hence, the main issue is to show that the integral 
is well defined on all of B2^ which is demonstrated here using the Caratheodory 
extension theorem and Dynkin's 7rA-theorem. 

Proof of Existence and Uniqueness for Theorem 3.1. Fix 1? G I. First we verify 
existence of a fluid model solution with initial measure d. For this, fix 1 < k < K 
and < t < 00. Given B £ 1Z, we have that B = [a, b) x [c, d) for some a, c £ R+, 
a < b < 00, and c < d < 00. Define / : [0, t] — > [0, 1] by 



m = (a+ w x r fe ) (B t _.). 



By (17), the fact that wi > and monotonicity properties of w(-), w(s) > for 
all s > 0. Then, since t > 0, we have that for s <= (0, t], 

f(s) = {s W ( S ) x r fc ) (B t - S ). 

Therefore, for s £ [0,t], 




t— s) - G k {d + t — s), if a + t - s < w(s) <b + t — s, 
otherwise. 



We see that a + t — s < w(s) < b + t — s if and only if a + t < w(s) + s < b + t 
if and only if r(a + t) < s < rib + t), where r(co) = 00. Hence, for s £ [0, t], 



f(s) 



G k (c + t-s)- G k (d + t—s), if r(a + t) < s < r(6 + t), 
0, otherwise. 



Then, since Gfc(-) and r(-) are continuous, / is Borel measurable and /„* f(s)ds 
is well defined. Further, since r(-) is monotone increasing, 

t i-tAr(b+t) 

f(s)ds= / G k {c + t-s)-G k {d + t-s)ds. 

JtAr(a+t) 

For B £ K, define 

lk {i){B) = MBt) + A fc J (<5+ s) x r fc ) {B t _ s )ds. 

It is clear that 7fc (*)($) = 0. Further, -f k (t)(M.\) < ■d k (M 2 + ) + \ k t < 00. So, 
in order to verify that Jk(t) is a finite premeasure on 1Z, it suffices to verify 
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countable additivity. This amounts to demonstrating that the summation and 
integral can be interchanged, which is an immediate consequence of the mono- 
tone convergence theorem. Then, by the Caratheodory extension theorem, 7^ (t) 
extends to a finite Borel measure 7^ (t) on B2 ■ 

We must verify that J%(t) satisfies (26) for all B £ 82- For this, we use 
Dynkin's 7rA-theorem [4, Chapter 1 Theorem 3.3]. Let 

V = {[a, 00) x [c, 00) : < a, c < 00}. (35) 

This is a 7r-system since it is closed under intersection. Let 

C = \b £ B 2 : 7j(t)(S) = M B t) + Xk J Q ( S Ms) x r fe) (Bt-s)ds 

Then £ is a A-system since £ £, C is closed under countable unions (by the 
monotone convergence theorem), and A\B £ C whenever B,A £ £ and B a A. 
Further, 

PcKc£c8 2 . 

Then, since the cr-algebra generated by V is B2, it follows from Dynkin's irA- 
theorem that C = 82- 

Since 1 < k < K and t > were arbitrary, it follows that 7* : [0, 00) — > , 
which is given by 

7*(i) = (7r(t),-,7Sr(*)). 
is a fluid model solution. Uniqueness is immediate since any fluid model solution 
C such that C(0) = i9 satisfies Ck{t){B) = l* k {t){B) for all B £ K (see [6, Theorem 
1.14]). □ 

Proof of Property 1 for Theorem 3.1. Fix § £ I. Let w(-) be the unique work- 
load fluid model solution with w(0) = and £(■) be the unique fluid model 
solution with £(0) = $. Since £(0) = <3 and w(0) = w$, uWo) = iu(0). Fix 
t £ (0,oo) and e £ (0, «;(*)). Set S = [io(t) - £,00) x R + . Then B( 2e ,o) = 
[w(i) +£,00) x R + . We wish to show that C+(t)(B) > and £ + (i) (B( 2s ,o) ) = 0. 
Then, since £ € (0, w(i)) is arbitrary, it follows that u?f(t) = 

Suppose that C+(i) (-B(2 E ,o)) > 0. Then by (26), there exists s € [0,i] such 
that w{t) + e + t — s < w(s), i.e., w(t) + t + e < w(s) + s. But w(-) + i(-) is 
strictly increasing, so that w(t) + t + e > w(s) + s, which is a contradiction. 

Thus, C+(t) (B(2e,0) )= 0. 

Since w(-) + t(-) is continuous and strictly increasing, there exists si £ [0,t) 
such that w(t)+t — e < w(si)+si. Then, for all s £ (si, t], w(t)+t — e < w(s) + s. 
Let S2 = (t — e?max) + - Then S2 £ [0, t) and t — s < d max for all s £ (s 2 ,t]. Let 
s* = s\ V s 2 . Then s* < t and 

K ft , N K ft K ft 

J2 / (<c w x r fc) (Bt- S )ds >J2 Gk ^ - s ^ ds -J2 Gk{t ~ s ^ ds > °- 



Hence, by (26), (+{t)(B) > 0. 



□ 
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The next goal is to verify Properties 2 and 3 for Theorem 3.1. For this, 
we state and prove the following lemma. Property 2 follows immediately from 
Lemma 4.1 by letting e decrease to zero. To verify Property 3, we must show 
that fluid model solutions change very little over short time intervals. Note that 
there are three mechanisms that cause the measure-valued function to change: 
fluid arriving, fluid departing, and the measure evolving. Verifying that the fluid 
departs in a smooth way is the main issue that needs to be addressed. We will 
use the result in the following lemma to assist with this as well. 

Lemma 4.1. Given i)el, let £(•) be the unique fluid model solution such that 
C(0) = For every T, e > 0, there exists k > such that 

,2^ sup sup (k(t)(C£ x } ) < e. 

l<fc<-K te [0,T] x,yeR+ 

To ease the reader's efforts to follow the proof of Lemma 4.1 given below, we 
outline the basic strategy. Having a good understanding of this deterministic 
argument will assist the reader in following the stochastic generalization used 
to prove Lemma 5.7. The basic idea in this case is to use (26) on a given k 
enlargement of a corner set. Then (24) can be used to bound the contribution 
from mass present in the system at time zero. Next one determines the time 
interval during which mass must arrive in order to contribute to the vertical 
portion of the enlarged corner set. As the reader will see, the end points of this 
time interval can be expressed in terms of the function t(-). Then one obtains a 
bound on the mass that can be present in the horizontal portion of the enlarged 
corner set. This second part turns out to be fairly easy to bound, as the reader 
will see. Hence the main difficultly is to bound the amount of mass that falls in 
the vertical portion. This relies on bounding the length of the aforementioned 
time interval, which can be done by demonstrating that the function w(-) + t(-) 
increases sufficiently quickly. But w(-) + t(-) actually increases quiet slowly at 
times when the workload fluid model solution is near d max , i.e., possibly at small 
times. So one must wait a short amount of time (until time 6 in the proof) during 
which just a small amount of mass enters the system before using this strategy. 
Once that small amount of time has elapsed, (16) can be used to obtain the 
desired bound. The mathematical details are given next. 

Proof of Lemma 4.1. Fix T,e > 0. Set A+ = £f =1 A fc , S = e/(4A+), M = 
w u V w(S), and c — J2k=i PkGk(M). Note that c > since M < c? max . Further, 
by monotonicity properties of workload fluid model solutions, w(u) < M for all 
u>5. 

For * € [0,T], 1 < k < K, x,y e R+, and K > 0, by (26), 

(c?., s) ) = a(o) ((^))J + a* jf (a+ w x r fe ) ((q;,,))^) 

By (24) there exists kq such that for all < k < kq, 

max sup sup Cfe(0) (fey)) J < j- 



ds. 
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Hence, for f 6 [0, T], 1 < k < K , x, y £ R+, and < k < k , 



Ut) (C{., v) ) < \ + X k / (d+ w x T k ) (C^ y) ) 



ds. 



(36) 



We must show that there exists k* < kq such that for all t £ [0, T], 1 < k < K, 
x, y £ R+, and < k < k* , the integral term is bounded above by 3e/4. 
Fix t £ [0,T], 1 < k < K, and x,y £ R+. For s £ [0, T] and k > 0, let 

/i fc (s, k) = G fe ((y -«;)++ 1 - s) - G k (y + K + t-s). 

For < k < kq, the integrand in (36) at time s £ [0,t] is bounded above by 

0, if w(s) + s < (x - k) + + t, 

1, if (x - k)+ + t < w{s) + s < x + n + t, 
hk(s,n), if x + k + t < w(s) + s. 

Equivalently, for < n < kq, using continuity and monotonicity properties of 
w(-) + l(-), the integrand in (36) at time s £ [0,t] is bounded above by 



0, if 0<s<t{{x~k)+ +t), 

1, if t((x - k)+ + t) < s < t(x + k + t), 
hk(s, t), if t(x + k + t) < s < t. 



(37) 



This together with (36) yields that for < k < kq, 
Ck(t) (C£, 



■,v) 



< i + A, 



r(x+K+t)f\t r t 

ds + A fe 

T((x~n) ++t)/\t 



T(x + K+t)/\t 



hk{s, n)ds 



< - + A fc (t(x + k + t) A t - t((x - k) + +t) At) 
rv+K 



(»-«)+ 



Gk(u)du 



< - + A fc (t(x + K + t)At- t{{x - k)+ +t)At)+ 2A/c(38) 

For < k < kq, let 

A(k) = t(x + k + t) A t - t((x - k) + +t) At. 

Fix < k < k . If t((x - k) + +t) >t, A(k) = 0. Also, if x + k + 1 < w(0), then 
A(k) = 0. Henceforth, we assume that t((x — n) + +t) < t and x + n + 1 > w(0). 
If t(x + k + t) At < 5, A(k) < 5. Otherwise, t(x + k + t) At > 5, and there are 
two cases to consider. First consider the case where r((x — k) + +t) > 5. Then 
(x — k) + + 1 > w(0) since t((x — k) + + t) > 0. Hence, by (29) and (16), since 
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t((x - k)+ + t) > 5, 

2k = X + K + t— (x — K + t) 

> X + K + t- ((X- k) + + t) 

= w(t(x + K + t)) + t(x + K + t) 

-w{t{{x - K) + + t)) - t((x - «) + + t) 
= G k (w(u))du 

f,—l Jt((x—k) 



r((x-K) + +t) 
K fT(x+K+t) 

> Y,P k G k (M)du 

k= l Jr((x-K) ++t) 



> cA(k). 

In particular, A(k) < 2k/c. The other case to consider is r((x — k) + + t) < 5. 
Then 

A(k) < t(x + K + t) At = (t(x + k + t) A t - 5) + 5. 
Further, by (29) and monotoncity of w(-) + t(-), 

x-K + t < (x - k) + + t < w(t((x - k) + + t)) + t((x - k) + + t) 

< w(5) + 5 < w{t(x + k + t)) + t(x + K + t) = X + K + t. 

Then, by (29) and (16), 

2k > x + K + t — w(S) — 5 

= w(t(x + K + t)) + t(x + K + t) - w{5) - S 



J i l-T(X-t-K-f-Z) 

= G k (w(u))du 
> VpJ G k (M)du 



k=l 

c(t(x + K + t) -5) 



In particular, 



9k 

A(«) <—+S, (39) 
c 

which is the largest of the four upper bounds. 

By combining (38), (39), and the definition of S it follows that for < k < kq, 

(Cf*,,)) < \ + ^ + \ k 5 + 2\ k K< t - + 2 -^ + 2\ k K. 
Let < k* < k be such that for all < k < k* 

max 2X k { — hi) K < 

Kk<K 
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Then, for all < n < «*, Cfc(*) ( C ( x . y )) < e - Since 1 < k < K,t e [0,T], and 
x,y,€ K+ were chosen arbitrarily, the result follows. □ 

Next, Lemma 4.1 is used to prove continuity for Theorem 3.1. For this, we 
need to specify a metric on M2 that induces the topology of weak convergence. 
For efficiency, we will specify such a metric on Mj for i — 1,2. Given i = 1,2 
and £, C G Mj, let d[£, £'] denote the Prohorov distance between ( and 
Specifically, 

d[C, C'] = inf{e > : ((B) < ('(B e ) + e and £(B) < ((B e ) + e 
for all closed B G Bi}. 

Then the following is a natural metric on Mf . Given i = 1, 2 and £, £' G Mj , 
define 

d K [(X'}=™x d[C k ,C k ]- (40) 

1</C<iv 

In the following proof, we show that for each fluid model solution ((■), lim/^o dx[C(i+ 
h),((t)] = for all t G [0,oo). 

Proof of Theorem 3.1 Property 3. Let 1) 6 I and let £(•) be the unique fluid 
model solution such that £(0) = §. It suffices to prove continuity of on 
[0,oo) for each 1 < k < K . Fix 1 < k < K. We prove continuity of Cfc(') on 
[0,T] for each T > 0. Fix T, £ > 0. By Lemma 4.1, there exists k such that for 
all < h < k, 

sup ( k (s)(C h )<s. 

se[o,T] 

Let < s < i < T be such that t ~ s < min(K, e/ max(2, Afc)). Note that 
y>w(u) x ( R +) = 1 ^ or an u e [0,oo). Let B <E B 2 be closed and set h = 
t — s > 0. Note that Bh C -B 2?i . By subtracting (26) applied to Bh at time s 
from (26) applied to B at time i, we obtain 

( k (t)(B) = Ck(s)(B h )+\ k f (6+ (u) xr k )(B t _ u )du 

< Ck(s)(B 2h ) + X k h 

< Ck(s)(B E ) + s. 

Consider (26) applied to B at time s, to B 2h at time t, and to C h at time s. 
For u G [0, s], we see that (w,p) G B u implies that (w — u,p — u) G B, which 
implies that either (w - u,p - u) e B C\ C h C C h or (w - u,p - u) e B \ C h . 
If (w - u,p — u) G C\ then (w,p) G C£. If (w - u,p - u) G B \ C h , then 
(w — u — h, p — u — h) G B 2 ^ and so (w,p) G B^ L +h . Hence, 

Ck(s)(B) < ( k (t)(B 2h ) + Ck(s)(C h ) < Ck(t)(B e ) + e 
Then for all < s < t < T such that < t - s < min(K, ej max(2, A fe )), 

d[Cfc(t) s Cfc(*)] <e. 
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Hence, Cfc(') is continuous on [0, T]. Since T > was arbitrary, Cfc(') is continuous 
on [0, oo). Since 1 < k < K was arbitrary, £(•) is continuous on [0,oo). □ 

Now that each statement in Theorem 3.1 has been verified, we prove Theorem 
3.3. 

Proof of Theorem 3.3. First suppose that $ € I is an invariant state. We must 
show that •& £ J, i.e., we must show that for some w £ [wi,w u ], $k(B) — 0™{B) 
for all B £ $2 and 1 < k < K. Since $ is an invariant state, Theorem 3.1 Prop- 
erty 1 implies that the unique workload fluid model solution such that u>(0) = 
is constant, i.e., w(t) — w$ for all t £ [0,oo). Then, by the monotonicity prop- 
erties of workload fluid model solutions, u>i < wd < w u . Let w — w§. Fix 
1 < fc < K, B £ V, and t > w. Then B = [a, oo) x [c, oo) for some a,c £ M.+ . 
Further, a + t — w > and (w — a) + < t. Hence, by (26), we have 

MB) = A fc / (6+ x T fc ) (B t _ s )ds 
Jo 

= A fc / Gfc(c + t- s)ds 

J (a+t-w)M 

= Afc / Gk{c + u)du 

Jo 

- V£{B). 

Since B £ V was arbitrary, i)k and 9™ agree on V . Since dk and 9% agree on V, 
they agree on 1Z and therefore they agree on Bi- So dk = 9^ . Since 1 < k < K 
was arbirary, -d = 9 W £ J. 

Next, fix wi < w < w u . Then the unique fluid workload solution w(-) with 
w(0) = w is constant, i.e., w(t) = w for all t £ [0,oo). Set = 9 W for all 
t £ [0, 00). Then, in order to show that £(•) is a fluid model solution, it suffices 
to verify that £(•) satisfies (26). For this, it suffices to verify that for all B £ V, 
1 < k < K, and t > 0, 

= ff£{Bt) + Afc /' (S+ x r fc ) (B t _ s )d S . 
Jo 

For this fix B £ V, 1 < k < K , and t > 0. We have 

WO + Afc / (6+ x r fc ) (B t _ s )ds 
Jo 

= 6^(B t )+X k f G k (c + t-s)ds. 

J(a+t-w)+M 

Case 1: Suppose that a > w. Then (a + t — w) + A t = t and it follows that 

6t{B t ) + Afc f (6+ x Tfc) (B t _ s )ds = = 0%{B). 
Jo 
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Case 2: Suppose that a < w < a + t. Then (a + t — w) + A t = a + t — w and 
9f(B t ) + X k f (6+xT k )(B t ^ s )ds 



o 

+ A fc / G k (c + t - s)ds 

J a-\-t—w 
w—a 



= Afe / Gk(c + u)du 
Jo 

= 8t(B). 

Case 3: If a + t < w, then (a + 1 - w) + A t = and 



t 

%(B t ) + \ k I (5+ x r fc ) {B t _ s )ds 



n 



= 6f{[a + t,w) x [c + t,oo)) + X k I G k (c + t-s)d 



is 

o 

pw — a—t pt 

= Xk Gk(c + t + u)du + Xk / Gk(c + u)du 

Jo Jo 

r-W — a r-t 

= Afc / Gk(c + u)du + Afe / Gk{c + u)du 
Jt Jo 

= Xk / Gk(c + u)du 
Jo 

= WW. 

Hence W is an invariant state. □ 
5. Proof of Tightness 

In this section, we prove that the sequence {Z™(-)}„ s n is relatively compact. 
For this, we apply [5, Corollary 3.7.4]. In particular, it suffices to prove compact 
containment and an oscillation inequality. This is done in Lemmas 5.5 and 5.11 
below. Throughout, 



K K 1 

X + = V X k and g + = V — 
^ ^ Ik 



k=l fc=l 

5. 1 . Preliminaries 

For 1 < k < K and t € [0,oo), let X* k {t) = X k t. For 1 < k < K, let £>£(•) = 
AJ(-)r fc . Define A*(-) = (X{ (■),..., A^(-)) and ©*(•) = (BJ (■)..., 2?£(.)). By a 
functional law of large numbers for renewal processes, asn-> oo, 



E"(-)^A*(-). 



(41) 
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Further, as a special case of Lemma 5.1 stated below, the following functional 
law of large numbers for the deadline process holds: asn-> oo, 

(V n (-),{ X ,V n (-)))^(V*(-),( X ,V*(-))). (42) 

To state Lemma 5.1, we need to introduce some additional notation and asymp- 
totic assumptions. 

Asymptotic Assumptions (AA) For each n € N suppose that we have 
K independent sequences of strictly positive independent and identically dis- 
tributed random variables that are independent of the exogenous arrival pro- 
cess. For each n e N, denote the fcth sequence by {g% j}ieN, and for each n e N 
and 1 < k < K, denote the distribution of 1 by rjj. We assume that TJ? has a 
finite mean 1/7/?, but we do not necessarily assume that F]J is continuous. We 
further assume that for each 1 < k < K, 

lim sup(xl(M.oo),n?) =0, (43) 

M->-oo neN 

and, as n — > oo, 

v n k % r fe . (44) 

Note that (43) is a uniform integrability condition and, together with (44), it 
implies that for each 1 < k < K, as n -— !• oo, 

« 1/7*. 

For each neN, let r n = T™ X • • • X T n K and let V = T l x • • • x T K . For n e N, 
1 < fc < if, and t € [0, oo), define 

<=i 

For neN, let £«(•) = • • • , S&O) and £«(•) = (0?Q, ,£&(•))• So 

then, for each neN, <?"(•) e D ([0, oo), Mf ) . We refer to £"(•), n e N, as a 
deadline related process. 

Lemma 5.1. Suppose that (AA) holds. Then, as n — > oo, 

(G n (-),(x,Q n (-))) ^ (v*(-),(x,-D*(-))) . 

Lemma 5.1 holds by (41) and [10, Theorem 5.1]. To see this, note the fol- 
lowing. Respectively, the number of classes and deadline distributions play the 
same role for the deadline related processes defined here as the number of routes 
and service time distributions do for the load process defined in [10]. Then [10, 
Theorem 5.1] holds under [10, Assumption (A)], and the conditions in [10, As- 
sumption (A)] relevant to [10, Theorem 5.1] are [10, (4.8)-(4.13)]. Conditions 
[10, (4.8)- (4.9)] hold here because the deadlines are assumed to be strictly pos- 
itive and have finite means. Condition [10, (4.10)] corresponds to (41) above. 
Conditions [10, (4.11)-(4.13))] hold here due to (43) and (44) above. 



Jennings and Puha/ Overloaded FIFO Queues with General Abandonment 



30 



For each n £ N and 1 < k < K, define the fluid scaled increments as follows: 
for < s < t < oo, 

E%(8,t) = E%{t) - E%{s) and \%{s,t) = \%(t) - \%{s), 

and 

Gt(s,t) = g2(t)-GUs) and Vl(s,t)=V* k (t)-VUs). 
Then, by (41), given T, e, r) > 0, 

liminfPf max sup \E%(s,t) - Xt(s,t)\ < e ) > 1 - r). (45) 

n^oo \l<*:<^o<s<t<T / 

Similarly, the following corollary holds as a consequence of Lemma 5.1 and 
continuity of the measures T k , 1 < k < K. 

Corollary 5.2. Let T,e,r) > and < a < b < oo. Suppose that (AA) holds. 
Then 



n— foo 



liminfP max sup \g%(s,t)((a,b)) - V* k (s,t){{a,b))\ < e ) > 1 - rj. 

K l<k<K <s<t<T 



We also wish to consider versions of the deadline related process for which 
the deadlines decrease at rate one until they expire. We refer to such a process 
as a residual deadline process. To define these, let {gki}ieN, 1 < k < K, and 
n £ N, satisfy (AA). For each neN, I < k < K and t € [0, oo), let 



and set 



Z?(0)+i 



E k(t) v n (t\ 

ntt) = E and - 



=i 



Then, for each n £ N and 1 < k < K, K%(-) £ D ([0, oo), Mi). For n £ N, 
set "£"(•) = (^(-),...,^(0) and £"(•) = (ft?(-)> • • • >^SrO)- For each 
1 < k < K and t £ [0, oo), let TZ>%{t) £ Mi be the measure that is absolutely con- 
tinuous with respect to Lebesgue measure with density \k (Gk(-) — Gf-(- + 1)). 
SetK*(-) = (Kt(-),...,K* K (-)). 

Lemma 5.3. Suppose that (AA) holds. Then, as n oo, 

K n (-)^K*{-). 

In order to prove Lemma 5.3, we must first prove that {7£™(-)}neN is tight. 
This is done using the same general approach as that outlined for proving tight- 
ness of {Z"(-)}„ e N- Specifically, we consecutively execute each step for proving 
tightness of both processes in Sections 5.2 through 5.4 below. Then in Section 6, 
we complete the proof of Lemma 5.3 by uniquely characterizing the limit points. 
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The result in Lemma 5.3 will then be used to help characterize limit points of 
{Z n (-)} n£ f,i and thereby help to complete the proof of Theorem 3.2. 

There are two specific choices of {gjtilieN, 1 < k < K and n G N, that will 
be of particular interest in what follows. 

Special Case 1: For n G N, 1 < k < K, and i G N, let — dk,i- This choice 
clearly satisfies (AA). For n G N, 1 < k < K, i G N, and tE [0, oo), let 

Ofc,»(*) = (dk,i - (t - tk,z»(p)+i) J ■ 
Then, for n E N and t E [0, oo), set 

■**(*) = E C Si (t) and = -•**(*)• 
i=i 

By Lemma 5.3, asn-> oo, 

(46) 

SpeciaZ Case 2; For n € N, 1 < k < K, and i E N, let g^. = d k ,i + v^. This 
choice satisfies (AA) as well. For n E N, 1 < k < K, i E N, and t E [0, oo), let 

«m(*) = + <i " (* " **,^(0)-h) + ) + ■ 

Then, for n G N and t G [0, oo), set 

Kit) 1 

»=i 

By Lemma 5.3, asn-> oo, 

V n (-)=^*(0- (47) 



5.2. Compact Containment 

In this section, we demonstrate the compact containment properties needed to 
prove tightness of {R- n {-)}nen and of {Z"(-)} nS N. In both cases, we will utilize 
[12, Lemma 15.7.5]. The application turns out to be simpler for {TZ n (-)} n eti 
since these are measures in M^. 

Compact Containment in Let K m = [0, to] and let denote its 

complement. Then {if m } m a forms a sequence of sets that increases to R+. By 
[12, Lemma 15.7.5], K' C Mf is relatively compact if and only if there exists 
a positive constant z and a sequence of positive constants {6 m }meN tending to 
zero such that for each ( G K' 

max Cfc(K+) < z and max Cfc(AT^) < b m > V to G N. 
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Note that for 1 < k < K and m £ N, C&C^m) < {x> Cfe) l m - Hence it suffices to 
show that there exist positive constants I and w such that for each £ £ K' 

max < ^ and max (v,Cfc)<u>- (48) 

We will verify (48) to prove Lemma 5.4. 

Lemma 5.4. Suppose that (AA) holds. Let T,rj > 0. There exists a compact 
set K C Mf suc/i i/iai 

liminf F(H n (t) £ K /or a// 1 £ [0, Tl) > 1 — 77. 

n—toc 

Proof. Fix T,r) > 0. For all f < fc < K and i G [0, T], 

^(t)(R+)<4"(T) and ( X , < ( X , ^(T)> . (49) 

Define 

0? = {m^ E%(T) < 2A+t| and flj = {max^ ( X , 9 n k (T)) < 2.9+t} . 

Set f$ = A? n By (45) and Lemma 5.1, 

liminf P(fij?) > 1 - 77. (50) 

n— >oo 

The result follows by combining (48), (49) and (50). □ 

Compact Containment in Mf- For each m £ N, let i£T m = [0, m] x [0, m] 
and let if^ denote its complement. Then {-fTmJmeN forms a sequence of sets 
that increases to By [12, Lemma 15.7.5], K' C M2 is relatively compact if 
and only if there exists a positive constant z and a sequence of positive constants 
{b m }m£N tending to zero such that for each (eK' 

max Cfe(Ki) < z and max (iK)<i m , VmeN. (51) 

l<fc<iC + l<k<K 

We will verify (51) to prove Lemma 5.5. 

Lemma 5.5. Let T, r\ > 0. There exists a compact set K C M^" smc/i that 
liminf P(Z"(<) £ K /or all t £ [0, T]) > 1 - 77. 

n— >oo 

Proof. Fix T, 77 > 0. First we identify a sequence of sets {SIq }„gN for which we 
will show that for each n S N on OJ, Z n (t) remains in a particular relatively 
compact set for all t £ [0, T]. By (27), there exists a compact set Ko such that 

liminf P(Z"(0) € K ) > 1 - 
n—>oo 4 
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Since Ko is compact, there exists a positive constant zq and a sequence of 
positive constants {a m } mg N tending to zero such that 

K C < C G My : max Cfe(Kl) < z and max C*C^™) < a™, V to G N 

l<fc<K + l<fc<if 

For neN, let 



n? = <^ max Z?. l (0)(Ri) < i and max Z?(0)(K^) < a m , V to G N 
[i<fe<K + i<fc<K R v M m/ ~~ 

Then 

liminfPfO?) > 1 - -. 



For n G N, let 



fi? = <^ max £?(T) < 2A+T 



By (45), 

liminfP(fiJ) > 1- -. 

Note that, for all n G N, 

JC E£(T) 

sup T^ n (o < r(o) + - y y v kti . 

Then, by (27) and a functional strong law of large numbers, there exists a 
positive constant u> such that 

liminfpf sup W n (t) < w ) > 1 - -. 



For n£N, let 



te[o,T] 



Og = { sup W n (t) < u) 



Notice that for each n, to G N, 1 < fc < K, and i G [0, T], 

(y f>?(T)) 

P£(i)((m, oo)) < 2)JJ(T)((m, oo)) < U ' fcl J/ ■ (52) 

TO 

Set c = 2A+5+T. For each to G N, let c rn = c/m. Then (52) and Lemma 5.1 
together imply that 

liminfP| max sup VI (*)((m, oo)) < Cm, V m G N ) > 1 - ~. 

n^oo \l<fc<#" te [0,T] / 4 
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For neff, let 

Q? = I max sup P?(*)((m, oo)) < c m , V to G N I . 
|i<fc<i^ e[0J1 " J 

Finally, for each n G N set 

njj = n 05 n ^3 n n£. 

It follows that 

liminfP(^) > I-77. (53) 

n— ^00 

Next we identify the relatively compact set K'. For this, let z = Zq + 2A+T 
and for to G N, let 



Z, 1 < TO < 10, 

am + Cm-wi m>w. 



Then {6 m } m eN is a sequence of postive numbers such that b m — > as to — > 00. 
Let 

K'= (6 Mf : max ( k (Rl) < z and max ( k (K c m ) < 6 m ,Vm G N I . 
[ i<fc<K + i<fc<x J 

Then, by (51), K' is relatively compact. 

It suffices to show that for each n G N, on Qft, Z"(t) G K' for all t G [0,T]. 
Fix n e N. On O™ fl 0£ we have that, for all 1 < k < K and t G [0, T], 

Z^)(]4) < Z£(0)(K 2 + ) + £»(T) < zo + 2A+T = z. (54) 

Next, by (13) under fluid scaling and the fact that for any x G M + and to G N, 
{Km)* ^ ^m- we have > for each m G N, 1 < fe < K, and t G [0, T], 



< - E ^Ki.pZj) 



n 
i=i 



A n k {t) 



j=Z»(0)+l 



Recall that on fi™, maxi<fc<^ Z^(0)(i4T^ l ) < a m for all to G N. Further, on fig , 
for to > t&, lr w „ >m | = for each 1 < fc < if and Z£(0) + 1 < j < A£(T). 

Additionally, on fig, for each 1 < fc < if and Z£(0) + 1 < j < A£(T), 
p2j < 4,i-z»(o) + W n {tl d ) < d fcJ -z»( ) + iZ>. 
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Hence, on n fij, for m > w and 1 < k < K, 

1 Am 

- E V >m) ^ ^fc(*)((»" - oo)) < c ™-- 

j=Zjf(0)+X 

Then on n fij n !!? for ra > ii, 

^(t)TO<a m + c m _«. (55) 

By (54) and (55), it follows that on Q^,_Z^(t)(K^) < b m for all m e N. This 
together with (54) implies that on fift, £ n (i) G K' for each i E [0, T]. □ 



5.3. Asymptotic Regularity 

This section contains results that are preparatory for proving the oscillation 
bounds. For the sequences of fluid scaled residual deadline processes and state 
descriptors, the sudden arrival or departure of a large amount of mass may 
result in a large oscillation. The focus here is showing that it is very unlikely 
that large oscillations take place due to departing mass. 

Asymptotic Regularity in Consider the sequence {TZ n (-)} n ^ of fluid 

scaled residual deadline processes. Given x € M+ and k > 0, let 

I£ = ((as - k) + ,x + k) . 

Note that for t € [0, oo), x £ R+, and k > 0, the mass in I£ at time t will all 
depart the system during the time interval (t + (x — k) + , t + x + n). This time 
interval is small if k is small. Hence, in order to avoid an abrupt departure of a 
large amount of mass, one needs to show that asymptotically such sets contain 
arbitrarily small mass. This is stated precisely in the following lemma. 

Lemma 5.6. Suppose that (AA) holds. Let T,e,r) > 0. Then there exists k > 
such that 



liminf P max sup sup ft£(i)(J£) < e) > 1 - 77. 

n ^°° \l<k<K t £[Q,T]xGR + J 

Proof. Fix T,e,rj > 0. Set k = §|- and M t = for < t < T. For each 

integer m > 3, n. € N, 1 < A; < K, and i € [0, T], 

0£(t)(((m - 2)«, 00)) < 0£(T)(((m - 2)«, 00)) < ^f^ - 

(m — 2)K 

It follows by Lemma 5.1 that for some m > 3, 

liminfpf max sup Ql (*)(((m - 2)k, 00)) < e ) > 1 - -. 

n-*-oo I l<fe<^ tG [o,T] / 2 
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Fix such an to and for n G N, let 



n? = \ max sup 0j?(i)(((ro - 2)«,oo)) < £ 

I l<k<K t g jq,T] 



Also, for neN, define 



Clo = i max sup max \Q2 (s,t)((jK, (j + 4)k)1 

\l<fe<if < s <f< T 0<j<m+M T 1 fcV AV ' V 7 " 

-^( a ,t)(0'K,0'+4)«))|< 



2M T 
By Corollary 5.2, 

liminfP(fig) > 1- -. 
For n£N, set fift = Q™ n fig. Then, we have that 

liminfPffi?) > 1 - r\. 

n—>oo 

Fix n G N, 1 < A; < K,x G R+, and < t < T. First, note that for any 
1 < i < El (t), we have ^ i+z » (0) < f, so that (t - t£ i+z „ (0) ) + = t - t% >i+Z n ioy 
It follows that 

Then for a; > wk and 1 < « < E%(t), 

(x + t — tj? j_|_^«(0) — K ) > a; — k > (m — l)/t > (to — 2)k. 



Hence, on fi™, for a; > tok 



< - E l((™-2)K,oo)(.9L) =0jf(t)(((m-2)K,oo)) <£. 



1=1 

Otherwise, a; < tok. Then 

Aft -2 BJ?(0'+i)«) 



Em 



i=BJ((Af t -l)«)+l fc .*+ z ?<°) 

Noting that for < j < M t - 1 and E%(Jk) < i < E%({j + 1)k Vt) we have that 

at + 1 - k - ifc, i+Z n (0) > a; + t - (j + 2)k > k 





X + t 




( 




-<->) 




K 
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and 



X + t + K - t% ii+Z n {0) < X + t - (j - 1)K < K ^ 



X + t 



■ 3 + 2 . 



Then it follows that 



M t -2 El((j+1) K ) 

nmw < - E E ^ 



j=0 i=E»(j«)+l 



(<?L) 



+- 



n 



E 



i /2 



i=E£((M t -l)K) + l 



(L^j-^-i,) 



(<?L) 



M t -2 



+CSf((M t -iKt) (^ m _ (Mt _ 1)} j. 

Since ir < ra/t, for all 1 < j < M t — 1, the left end point of /^ x+t j _^ satishes 



< K 



j - 2/t <i + i<mK + M t K < (m + M t )k. 



Also note that < t — (M t — 1)k < t — (t/n — 1)k = ft. Then, on £1%, since 
x < run, 



M t -2 



+2>i((ilf t -l)M) (^(L^j-^-D^ + l 



< 



M t -1 

E AfeK 



F k K 



x + 1 



.7 + 2 
x -\- 1 



< \ k K e [^(0' + 4 )«) - + 

J=-2 



3 -2 

£ 

2' 



So then, using the fact that ffc(-) is a cumulative distribution function and that 
each point in R + is included in at most four intervals of the form (jn, (j + 4)k], 
j = —2,-1,0,1,2, Thus, on QQ, since x < ran, 

<4A feK + | <£. 

Since, n e N, 1 < k < K , x e M + , and t e [0,T] were chosen arbitrarily, this 
concludes the proof. □ 
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Asymptotic Regularity in We need to prove an analog of Lemma 5.6 

for {Z n (-)} ne jq. In particular, we wish to prove the following prelimit version of 
Lemma 4.1. 

Lemma 5.7. Let T,e,rj> 0. Then there exists n > such that 

liminf P f max sup sup Z£(t)(C£) < e ) > 1 - 77. 

n ^°° y^ k ^ K xm%te[o,T] J 

Before proving Lemma 5.7, we verify the following regularity result for the 
initial state, which is the stochastic analog of (24). 

Lemma 5.8. Let e, r) > 0. Then there exists k > such that 

liminf P f max sup Z£(0)(C£) < e) > 1 — 77. 
n->oc \i<k<K xeR 2 + J 

Proof. Fix e, 77 > 0. Given i = 1,2, recall definition (25) of the projection 
mapping m : M2 — > Mi. We apply the argument given in [8, Pages 835-836] for 
measures in Mj to the projection mappings applied to Z™(0), n G N, to verify 
that there exists k > such that 

liminfpfmax sup (l [(ic _« ) +, :c+K] , tt, (Z™(0))) < J J > 1 - r,. (56) 

The desired result follows from (56) since for all n G N, 1 < k < K, x = 
(xi,x 2 ) G K+ and k > 0, 

^fc(°)(C?) < (iKxi-^ + .ai+re]' 71 "! (^+(°))) + ( 1 [(K 2 -K)+,a;2+«] I 7r 2 (-Z+(0))) . 

In order to verify (56), we must verify that suitable i^-dimensional analogs 
of [8, (3.19)- (3.22)] hold. For this, for ( G Mf and i = 1,2 we adopt the 
shorthand notation 

MO = (MCi), ■ • ■ , MCk)) and (xMO) = ((X, ^(Ci)> , . . ., (xMCk))) . 

Note that for v G M2 and i = 1,2, (x> 7r i( 1/ )) — (x°PiT l/ ) — {Pii v )- Then, by 
(27), as n — > 00, 

(^(^"(O)),^^"^)),^,^^"^))),^,^^"^)))) 

=>■ ( 7I "l(- Z o)»7 I "l(- Z o)><X> 7I "l(- Z o)) ><X,7Tl(^))), 

which is the if-dimensional analog of [8, (3.19)]. Using (A. 3) we obtain the 
following i^-dimensional analog of [8, (3.20)]: 

max maxE Ul^AZZ k ))] < E [Z * ,(R+)1 < oo. 

Using (A. 2) we obtain the following i^-dimensional analog of [8, (3.21)]: 
max max E [(x,^^*^))] <E[{ Pl + p 2 , Z* + )] < 00. 

Kfc<i\ 1 — 1,2 
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Using (A.l) (and in particular (1.1)) gives the following if-dimensional analog 
of [8, (3.22)]: 

P SU P (!{*}> 7Ti(£o,fc)) = =1. □ 

n<Kifi=i,2 l6R+ i 

Before moving on to prove Lemma 5.7, we obtain an almost sure upper bound 
on the mass in C* x y ^ in an arbitrary coordinate of the nth system at time t for 
each k > 0, x,y G E+, and t G [0, oo). To this end, let n G N, 1 < k < K, 
t e [0,oo), x,y £ R+, and k > 0. The nth system analog of (15) for the set 

Am 



j=Z»(0)+l 

But then, since the (a:, y)-shift of a set followed by the K-enlargement contains 
the K-enlargement followed by the (x, y)-shift, 

^W(c , rx, a) )<^(o)(c'f x+t)t/+t) )+ E («,;(*))• ( 5? ) 

j=Z»(0)+l 

We simplify the summation term by focusing on each coordinate separately. For 
this, we classify jobs by those whose residual deadline time causes the associated 
unit atom to lie in a certain horizontal band and those whose residual sojourn 
time causes the associated unit atom to lie in a certain vertical band. Then, for 
each n G N, 1 < k < K, t G [0, oo), x,y £ R + , and k > 0, 

zm(C^ y) ) < Z[!(0)(C^ y+t) )+ J2 li.-Kj(*)) (58) 

j=2»(0)+l 

Afe(t) 

+ E 

We are prepared to prove Lemma 5.7. The proof given below can be regarded 
as a stochastic version of the proof of Lemma 4.1. 

Proof. Fix e, rj > and T > ^j^- We begin by defining a sequence {f2Q} rag N 
of events on which the prelimit processes satisfy properties analogous to those 
exhibited by fluid model solutions and used in the proof of Lemma 4.1. By 
Lemma 5.8, there exists kq > such that 

liminfP [max, sup Z fe "(0) (c£ y) ) < - 1 > 1 ' 



n— foo 



A ~ 6 
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For neff, let 

fi? = J max sup J?(0) (c?° ,) < 7 I . 
By Lemma 5.6 (see special cases 1 and 2), there exists a k\ > such that 
liminf P ( max sup sup A n k {t){I^) < - j > 1 - | 



and 



liminfPf max sup sup V^ftYH 1 ) < - ) > 1 - -. 

n-KM \ l<fc<^ aG R +t6 [ ,T] 4/ 6 



For neN, let 



^ = max sup sup A n k (t)(I^) < j 



0?=J max sup sup V fc n (t)(L fil )<^ 



and 



For neN, let 



flj = max sup A2(t)-A2(»)<2A+(t-«) 

^l<fc<-K o<s<t<T 

Since for all n G N, 1 < fc < K, and < s < t < T, A%(t) - A k l (s) = 
E%(t) - E%(s), (45) implies that 

liminfPf) > 1- -■ 

Next we identify positive constants 5 and M, analogous to the constants 6 
and M defined in the proof of Lemma 4.1. Let S — 24 ^ + . Then S < T. In 
order to define M, there are two cases to consider, based on the nature of the 
abandonment distributions. 

Case 1: First suppose that d max < 00. Let w max (-) denote the maximal workload 
fluid model solution, i.e., the workload fluid model solution such that u> max (0) = 
d mH . By the relative ordering property of workload fluid model solutions, it 
follows that for any workload fluid model solution w(-), w(t) < w max (i) for all 
t G [0, 00). In particular, for any workload fluid model solution w(-) for all t> 8, 

w(t) < 

^max 

(S). (59) 
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Set M\ = (w max (5) + d max ) /2. By monotonicity properties of workload fluid 
model solutions, w u < u> max (<5) < Mi < d max . It follows from (27), (A.l), 
Theorem 2.2, and (59) that 



P sup W*[t) < Mi = 1. 
\t&[6,T] J 

Case 2: Next consider the case d max = oo. By (27) and (A.l), there exists 
Mi > w u such that 

F(W*(0) < M 2 ) > 1- -. 

6 

Since M 2 > w u , (27), (A.l), Theorem 2.2 and monotonicity properties of work- 
load fluid model solutions imply that 

P ( sup W*(t) < M 2 ) = P (W*(0) < M 2 ) > 1 - \. 
\te[o,T] J 6 

Set 

M 



Mi, if<i max <oo, 
M 2 , ifd max = oo. 



Now we proceed to bound the prelimt processes by M with probability asymp- 
totically close to one. For n € N, let 



ft™ = i sup W n {t) < M } . 
[te[8,T] J 

Note that by (27), (A.l), Theorem 2.2 and that fact that workload fluid model 
solutions are continuous, the convergence in distribution in (21) takes place with 
respect to the topology of uniform convergence on compact sets. Furthermore, 
the set {/ £ D([0, oo),R+) : sup tg j,5 T ] f(i) < M} is open with respect to this 
topology. Then by the Portmanteau theorem 

liminf P(fi?) > 1 - -. 

Set 

k , e ec 

c= / PkGk(M) and k = min I kq, ki, 



^k~*^j — V 24A+'72A H 

Note that c > since M < d max , and so K > 0. Finally, for n € N, let 

"6 = 1 ™p |X"(t)|<^l. 

^te[o,T] z J 

Then, by (20), 

liminfP(fi") > 1- -. 

n— >oo 6 
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For neN, set 

n 7 > = n n£ n Q£ n o£ n n nj. 

Then 

liminfP(fi£) > I-77. (60) 

Fix n € N such that 1/n < e/12, 1 < fc < K, x, y € K+, and f e [0,T]. The 
fluid scaled analog of (58) is 

« ^(*) 

(<?.„,>) < 4 n (o)(^ +t , a+f) ) + - £ ^K,W) (si) 

j=Zjf (0)+l 

+- E i^(pSj(*))- 

j=Z»(0)+1 

We will show that on J7q , the right hand side of (61) is less than or equal to e. 
Since k < «o, we have C£ x+t y+t) C C£ +tjl , +t) . Then, on Of, 

4"(0) (Cf x+ti3/+t) ) < |. (62) 

Also, for all Z£(Q) + l<j< A%(i), the quantity p% d (t) is either a^._ z „ (0) (t) 
or v k j-z n (o)(t)- ^ follows that, on £1% n O3 , 

j=Z™ (0)4-1 

Combining (61), (62), and (63), we see that on £7q , 

3?(t) < ^ + - E Ifl «*(*)) • ( 64 ) 

j=ZjJ (0)4-1 

Finally, we bound the second term on the right side of (64). To this end, 
define the prelimit version t™(-) of r(-) as follows: 

r n (s) = inf {u > : W n {u)+u > s} , s£ [0,oo). 

Notice that for all s E [0, 00), W n (s) + s > s, so that t"(s) € [0,s]. During 
busy periods, W n (-) + *.(•) jumps up exactly when jobs arrive that contribute 
to the workload, and remains constant otherwise. During idle periods, W n (-) + 
(,(•) increases at rate one. In particular, W n (-) + /,(•) is right continuous and 
nondecreasing. Then, 

W n (T n (s)) + T n (s) > s, for se [0,oo), (65) 
W n (T n (s)-) + T n (s) < s, for s e (W n (0),oo). (66) 
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The function t™(-) will be used to determine which arrivals may contribute to 
the sum in (64). First we show that 

, 

~ E ^jKiW) < A n k {T n {x + K + t)M) (67) 

-Al (r™((z -«)++*) At) +-. 
To see this, let Z™(0) + 1 < j < A%(t), Then t% d < t and 

(x - k) + < W%j (t) < x + k +> (x - k) + + t < W k \j + t^ j < x + K + t. 

But w% d = W n (tl tj ). Then, since W n {-) + l(-) is nondecreasing, t% tj < t and 
w k,j (*) e I x imply that 

r" Ka: - k) + +() A(< tjjj < r"(x + k + t) A t. 

Therefore, (67) holds. 

Next, we proceed to bound the right hand side of (67) on 17q. By (64), (67), 
the definition of f2J, and the fact that n is such that 1/n < e/12, it suffices to 
show that on fig , 

2A+ (t' 1 {x + n + t) At- T n ((x - k)+ + t) A t) < %. (68) 

6 

If T n ((x — k) + + t) > t, then the left side of (68) is zero, and so (68) holds. 
Henceforth, we assume that r n ((x — k) + +t) < t. If r n (x + n + t) A t < S, then 
(68) holds since 5 = e/24A + , which implies that the left hand side of (68) is no 
larger than e/12. Otherwise, 6 < T n (x + n + t)At (so that r n (x + K + t) > and 
then x + K + t > W n (0)). First consider the case where S < r n ((x — k) + + t). 
Then, by (65), (66), and the nondecreasing nature of W n (-) + t(-), 

2k = x + k + t — [x — K + t)>x + K + t— ((x — k) + + t) 

> W n (T n (x + K + t)-)+T n (x + K + t) 
-W n {T n ((x - K)+ + t)) - T n ((x - K) + + t) 

> W ,l (T n (x + K + t) A t-) + T H (X + K + t) At 

-W n (T n ({x - k)+ + t)) - T n ((x - K ) + + t). 
Using (19) and the nondecreasing nature of the idle time process, we obtain 
2k > X n (T n (x + K + t)At-)-X n (T n ((x-K)+ + tj) 

K pT n (x + K + t)At 

+ Y,P k G k {W n (u))du. 

k=l Jr"-((x- K ) + +t) 

Since 8 < T n ((x - k)+ + t), on n 

2k>-k + c (r n {x + k + t) A t - T n ((x - k) + + t)) . 
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By definition of k, 3k/c < e/24A+, which implies that the left side of (68) is no 
larger than e/12 on fig n O5 , and so (68) holds. The last case to consider is the 
case T n ((x - k) + +t) < 5 < r n (x + K + 1) A t. Then, by definition of 6, on fij, 

2A+ (r n (x + K + t)At-T n ((x- n)+ + 1) A t) (69) 

~ T2 + 2A+ + K + t ') At - S )- 

By (65), (66), and monotonicity of W n {-) + <,(•), 

x — K + t < (a; — k) + + t 

< W™(r n ((a; - k)+ + t)) + r n ((x - k)+ + t) 

< W n (S) + S 

< W n (r n {x + K + t) A t-) + T n (x + K + t) A t < x + n + t. 

This together with (19) and the nondecreasing nature of the idle time process 
implies that, on Of n Hg, 

2k > W n (T n (x + K + t)At-)+T n (x + K + t)At-W n (5)-6 

> X n (T n (x + n + t) At-) - X n (6) 

+ Vp fe / G k (W n (u))du 

> -K + c(T n (x + K + t) At- 5). 

Using the definition of n and (69) implies (68). □ 



5-4- Oscillation Bounds 

This section establishes the second main ingredient for proving tightness of 
the fluid scaled residual deadline processes and state descriptors, controlled 
oscillations (see Lemmas 5.10 and 5.11). Both proofs are similar in spirit. The 
one for the residual deadline processes is slightly simpler, so it is presented first. 
For this, recall the definition of dx given in (40). Then the modulus of continuity 
is defined as follows. 

Definition 5.9. Let i = 1,2. For each £(•) e D([0, 00), Mf ) and each T > 5 > 
0, define the modulus of continuity on [0, T] by 

w T (C(-),<5)= sup sup d K [C(t + h)X{t)]. 
te[o,T-S] he[o,s] 

Lemma 5.10. Suppose that (AA) holds. For all T > and e,n € (0,1), there 
exists S € (0, T) such that 

liminfP(w T ("K™(-),<5) < e) > 1 - r). 
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Proof. Fix T > and e,r) £ (0,1)- By (41) and Lemma 5.6, there exists a 
K £ (0, e) such that for any 5 £ (0, T), the events 

fi? = I max sup fcjj(*)([0,«])<4, 

ft£ = I max sup sup + ft) - Ejf (t)] < 2A + 5 1 , 

[i<fc<-K ie[o,T-5] he[o,«] J 

satisfy 

liminfP(Q£) > 1 -n. (70) 

n— >c~ 

Fix such a k and set S = min ( k, 



2A+ 

Fix n G N, 1 < k < K, and a closed set B £ B x . Let t £ [0, T - 5], ft e (0, 5], 
and 1 < i < E%(t). If g^(t) > ft or ^(t + ft) > 0, then 

9k,iV)=9k,i(t + h) + h. 
Then, for ft £ (0, <5], since h < 6 < k < e, 

n(t)(B) < K n k (t+h)m+n n (t)({o, K }) 

K n k (t + h)(B) < ni{t)m + El{t + h)- E%(t). 

Then, on Qft, for all i 6 [0, T — £] and ft G [0, <5], 

d(ft£(t + ft), <e. 

This together with the fact that n£N. 1 < k < K, and the closed set B e 6 2 
were arbitrary and (70) implies the result. □ 

Next we generalize the preceding argument to prove the following analogous 
lemma for the fluid scaled state descriptors. There are two main distinctions. 
One is that k < e/2 since in ft > time units a unit atom moves along a diagonal 
path a distance of \/2h. The other is that (15) has been established. 

Lemma 5.11. For all T > and e, n G (0, 1), there exists 5 G (0,T) such that 
liminf P(w T (Z™(-),<5) < e) > 1 — 77. 

n— >-oo 

Proof. Fix T > and e, rj £ (0, 1). For each k > 0, let C K be the closure of C K , 
or 

(7 K = [o, K ] x M+ UM+ x [0,k]. 
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By (41) and Lemma 5.7, there exists a k £ (0, e/2) such that for any <5 £ (0, T), 
the events 

fi? = I max sup 4"(t)(C K )<el, 

^l<fc<if tg[ ,T] J 

fij = | max sup sup [££(i + /i) -££(*)] <2A+<5 1, 
|^i<fe<K t e [o,T-$] fce[o,«] J 

^ = n^nnj, 

satisfy 

liminfP(O^) > 1 -r). (71) 

n— >oo 

Fix such a k and set 5 = min (k, jx^J ■ 

We begin by noting two basic facts that will be used in the proof. Firstly, for 
all B e B 2 and h £ [0,(5], 

BC{B e ) h UC K . (72) 

To see this, take some B £ B 2 , h £ [0,6], and (w,p) £ B\C K . By the construction 
of C K , we have w,p > k > <5 > h. Because h < S < n < e/2, it follows that 
(w — h,p — h) £ B e and (w,p) £ {B £ )h- In addition, since 6 < e/2, we have that 
for all B £ B 2 and h £ [0,(5], 

B h C B c . (73) 

Fix n £ N, 1 < k < K, and a closed set B £ B 2 . Let B = B \ C. Note that 
£ € 6 2 ,o and by (11), Z%(t)(B) = Z%{t){B) for all t £ [0,T]. Then, we can 
then conclude from (11), (15), and (72) that, on fift, for all t £ [0,T - 5} and 
he [0,(5], 

^?(t)(B) = mm 

< Z%{t){{B%) + Z%{t){C«) 

< Z%{t + h){B*) + e 

< Z%{t + h){B*)+e. (74) 

Also, by (11), (15), (73), and the fact that 5 < ^r— , it is true that on Qq, f° r 
allte [0,T-(S] and /i € [0,(5], 

4"(i + />)(£) = + 

< zmiB h )+m+h)-m) 

< Z%(t)(B h )+e 

< Zmm+e. (75) 
Because k and -B were chosen arbitrarily, (74) and (75) imply that on SIq 

w T (Z n (.),6)<e, 

The result follows from this and (71). □ 
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6. Characterization of Limit Points 

The main goal of this section is to prove Theorem 3.2. First we prove Lemma 
5.3, which is then used in the proof of Theorem 3.2. 

Proof of Lemma 5.3. Throughout, we assume that (AA) holds. Together Lem- 
mas 5.4 and 5.10 imply tightness of {72™(-)}neN. Let M C N be a strictly increas- 
ing subsequence tending to infinity and 7£(-) a process such that as m — > oo, 

n m {-) =►£(•). 

By Lemma 5.6, lt(t) doesn't charge points for all t £ [0, oo) almost surely. By 
(41) and Lemma 5.1 and the deterministic nature of those limiting processes, 
asm^oo 

(^™(.),S m (.),a m (-),<X^ m (-)))^(^(-),^(-),^(-),(x,^(-)>)- (76) 

Using the Skorohod representation we may assume without loss of generally that 
all random elements are defined on a common probability space (fi, J 7 , P) such 
that the joint convergence in (76) is almost sure. Fix lu £ il such that TZ(t)(u) 
doesn't charge points for all t £ [0, oo) and as m — > oo 

(£ m (-)M, £ m (-)M, £ m (-)M, <x, G m (-)) M) (77) 

Henceforth, all random variables are evaluated at this lu. It suffices to show that 
U(-) = K*(-). For this, it suffices to show that for all 1 < k < K, t £ [0, oo) and 

x £ R+, 

H k {t){x,oo) = A fc f G k (x + t-s)ds. (78) 
Jo 

To see this, note that for all 1 < k < K, t £ [0, oo), and x £ K+ 

/•OO 

n% (t) (x, oo) = A fc / [G k (y) -G k (y + t)} dy 

J X 

rX-\-t rt 

= Afc / G k (y)dy = \ k / G k {x + t- s)ds. 

Jx JO 

Next we verify (78). For this, fixed 1 < k < K , t £ [0, oo), and x £ R+. Given 
L £ N, let k = t/L and set 

t( = Ik and xg, = x + 1 — ti, for i = 0, . . . , L. 

Then, given L £ N, 1 < i < E%(t) if and only if there exists I £ {0, . . . , L - 1} 
such that t™ Z m( ) +i £ (t£,t e+ i}. Further, given L £ N, 1 < i < E%(t) and 
t £ {0, ...,L — 1} such that i™ Z m( ) +i € (i^>£e+i], < 5™^ implies that x < 
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g%i(t)- Similarly, given L G N, 1 < % < E%{t) and I G {0, — 1} such that 

*™z™(o)+i G (*«>*m]) x < implies that x e+1 < gj™ . Hence, given L G N, 



L-l 



L-l 



By (77), given LeN and e > there exists M G N such that for all m> M 
max (tt,te + i) (xt,oo) - V* k {tt,tt+\) {xt,oo)\ < 



o<e<L-i 



2L' 



max \C™ (t e ,t e+1 ) (x£ + i,(X)) -V* k (t e ,t£ + i) (x e+ i,(X))\ < — 



0<£<L-1 



2L 



Hence, given LgN and e > there exists MeN such that for all m > M, 

L-l L-l 

£ -<KT(t)(x^)<J2- K 2 



X^^fc (t e ,t e+1 )(x e ,(X)) - ^ <Tl™(t)(x,00) < ^2 V k (U,U+l) (x i+ i,00) + £ 



t=0 

Given L G N, we have that 

L-l 



1=0 



L-l 



X^^fc (tt,t e+1 ) (x t ,oo) 

1=0 
L-l 



L-l 



Respectively these are upper and lower Riemann sums, and since is con- 

tinuous, they both converge to J* Gfc(x + t — s)ds as L —> oo. Given e > 0, 
let L G N be such that 



L-l 



L-l 



f=0 



£=0 



e 

<2' 



Hence, given e > 0, there exists M G N such that for all m > A/, 



R%(t)(x,oo) - A fe / G fc (a; + i- s)ds 



< e. 



□ 



Thus, (78) holds, as desired. 

Having proved Lemma 5.3, we are now ready to prove Theorem 3.2. 

Proof of Theorem 3.2. Together Lemmas 5.5 and 5.11 imply tightness of {Z n (-)} n ^. 
Let M C N be a strictly increasing subsequence tending to infinity and Z*(-) a 
process such that as m — > oo, 



Z m {-) => Z*{-). 



(79) 
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Note that by (A.l), Z*(-) £ I almost surely. Hence, in order to prove Theorem 
3.2, it suffices to show that Z*(-) satisfies (26) almost surely. Indeed, once this 
is verified, it follows by the uniqueness asserted in Theorem 3.1 that the law of 
the limit point Z*{-) is unique, and so Z n {-) =$> Z*(-) as n — > oo. 
By (27), (A.l), and Theorem 2.2, as m -> oo, 

W m {-) => W*(-), (80) 

where W*(-) is almost surely a workload fluid model solution such that W*(0) 
is equal in distribution to Wq. We would like to argue that this convergence is 
joint with (79). By (27), (20), (46), (47), and the fact that the limit in (20) and 
1Z*(-) are deterministic, as m — > oo 

(z m (-),w m (o),x m (-),A m (-),v m (-)) (z*(.),w ,o,n*(-),n*(.)). 

This together with (22), (23), and (80) implies that as m — > oo, 

(Z m (-),W m (-),X m (-),A m (-),V m (-)) (81) 
=4- (Z* (■),]¥*(■), 0,11* (■), 11* (■)). 

Using the Skorohod representation we may assume without loss of generally 
that all random elements are defined on a common probability space (fi, T , P) 
such that the joint convergence in (81) is almost sure. By (27), Z*(0) satisfies 
w z-(o) — W*(0), (A.l), (A. 2), and (A. 3) almost surely. Furthermore, by Lemma 
5.11, Z*{-) is continuous almost surely. In addition, by Lemma 5.7, 

P (Z* + {t){C x ) = for all t € [0, oo) and x £ R"j_) = 1. (82) 

(cf. [9, Lemma 6.2]). Fix lu E fl such that W*(-)(ui) is a workload fluid model 
solution and Z*(-)(ui) is continuous and satisfies (82), tUz*(o)(w) — W*(0)(u)), 
(A.l), (A. 2), (A. 3) and, as m -> oo, 

(z™(.)M,iu m (-)H^ m (-)M,v m (-)M) 

^(Z*(.)(w),W*(.)(w),1l*(.),1l*(.)). 

For all t e [0, oo) and m £ M, set 

C m (t) = Z m (t){u) and u; m (t) = W^ m (t)(w), 

({t) = Z*{t){uj) and = W*(t)(u). 

By (27), w(0) = w$ where $ — C(0). In order to prove Theorem 3.2, it suffices 
to show that £(•) is a fluid model solution for the supercritical data (A, fx, T) and 
initial measure fl. In particular, we must show that £(•) satisfies (26). In this 
regard, recall that V is a 7r-system (see (35)). As in the proof of Theorem 3.1 
given in Section 4, it is enough to show that £(') satisfies (26) for all B £ V . 

Fix B £ V . Then B = [a. oo) x [c, oo) for some < a, c < oo. Fix t £ [0, oo). 
In what follows, all random elements are evaluated at the specific uj fixed in the 
preceding paragraph. Since £(•) is continuous, C™( s ) — ► C( s ) f° r an s G [0,^]- By 
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(82), C+(0)(C (a+t , c+t) ) = and ( + {t)(C (a , c) ) = 0. Then, for each 1 < k < K, 
we have 

lira Ck (0) (B t ) = C k (0) (B t ) and lira Ck (*) ( B ) = ( k (t) (B) . 

However, by the fluid scaled version of (15) with h = t and t = 0, for each 
m G M and 1 < k < K, 



1 ^r(t) 



Hence, in order to verify that (26) holds, it suffices to show that for each 1 < 
k < K, 



lim 

m^-oo 777, 



Note that for 1 < fc < K, 



J (fw(a) X Th ) ( B t-s)ds = Afe. y l{ w ( s )> a +t-s}Gk{c + t - s)d,S 



= Afe / Gk(c + t~s)ds 

Jr(a+t)At 
pc+t—r(a+t)At 

= Afe / G k (u)du. 



Hence, in order to verify that (26) holds, it suffices to show that for each 1 < 
k < K, 

lim - f t lB t ^(w m ( S ),pT, Ar{s) )dEr(s) (83) 

m->oo 171 Jq V fc / 

/■c+t— T(o+t)At 

= Afe / G k (u)du. 



In order to verify (83), fix 1 < k < K. Note that for all m G M and 1 < i < 

d k ,i < Pk,zr<to)+i ^ <i- 

Then 

i Bt _. (^ m ( S ),4 )Brw )^r(s) 

U_. ( w m ( s ),^ r(s) )d£;fe m ( s ) 

(>(*), 4,^) +<^ (s) ) dSTW- 



< 



< 
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Let < 6 < t. Take M such that for all m > M, sup < s < t \w m (s) - w(s)\ < 5. 
Then, for all m > M, 



< 



< 



ls t _ 8 (w{s)-5,d ktET(s) )dE%{s) 

l Bt _ s (w m (s),pZ ATis) )dE?( S ) 

l flt _. (w(s) + S, d Ws) + < £r(s) ) dE?(a). 
For s £ [0, oo), w(s) — S > a + t — s if and only if s > r(a + t + 5). Then 

(ttf(*)-J,d fc , B m W )£iB^(*) 

(d fc|Bjrw ) dsr(s) 

l[c,oo) (4,B-( S ) - (* - S)J cLE^O) 
L[c,oo) (< £r(s) W)^ m (s) 



(a+t+5)At 
t 



1 [c+£ — s, oo) 



T(a+i+<5)At 
t 

lr 



r(a+t+S)At 

lr. 



[coo) (o^w^iETW 
= ^(*)([c, oo)) - A^(r{a + t + 6) A *)([c + i - r(a + * + <5) A i, oo)). 

Similarly, 

ls t _ s (w(s) + ,5, d fc , B m W + »^ B m (a) ) dE^(s) 
= V£*(t)([c, oo)) - Vf(r(a + i — <5) A *)([c + 1 - r(a + 1 - 8) At, oo)). 
Then, by Lemma 5.3, 

i r* 



lim inf 

m— too ui 



> 7ifc(t)([c, oo)) - 7e£(r(a + f + S) A <)([c + t - r(a + i + (5) A i, oo)) 

But 

7^(*)([c, oo)) - 7^(r(a + t + S) A t)([c + t - r(a + t + 5) A i, oo)) 
= A fc / (G k (u) - G k (u + t))du 



-X k / (G k (u)-G k {u + T (a + t)At))du 

J c+t-r(a+t+S)At 
i-c+t-r(a+t+S)At 

X k / G k (u)du. 
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So then, since 5 £ (0,i) is arbitrary, letting S \ and using continuity of r(-) 
yields that 



1 /•* / . fc+t— r(o+t)At 

liminf- / l Bt _ s U m (*),P^ r(s) dBT(«) > A fc / G fe («)cfa 



/•* z x fc+t— r(o+t)At 

m 

Similarly, 

limsup- / l Bt _, ( W m (s),^ Ar(5) dS fc m (s) < A fe / 



Gk{u)du. 



Hence (83) holds, as desired. □ 
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